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EXECUTIVE  SUMMRY 


This  document  serves  as  a  user's  guide  for  the  Multiple  Aircraft 
Instantaneous  Line  Source  (MAILS)  atmospheric  dispersion  model  computer 
software.  MAILS  is  an  air  quality  screening  model  that  provides  conservative 
estimates  of  ground-level  pollutant  concentrations  resulting  from  aircraft 
engine  emissions  along  low-altitude  (under  3000  ft)  military  training  routes 
(MTRs). 

The  MAILS  model  was  developed  to  evaluate  the  air  quality  impacts  of 
proposed  MTRs,  in  accordance  with  the  environmental  impact  analysis 
requirements  of  the  National  Environmental  Policy  Act  of  1969  (NEPA,  Public 
Law  91-190)  and  AFR  13-2,  USAF  Environmental  Impact  Analysis  Process. 

Previous  analyses  indicated  that  the  air  quality  impacts  of  low-flying 
aircraft  were  potentially  significant  with  respect  to  Prevention  of 
Significant  Deterioration  (PSD)  Class  I  pollutant  increments  (applicable 
primarily  in  certain  national  parks  and  wilderness  areas)  and  were 
insignificant  with  respect  to  National  Ambient  Air  Quality  Standards  and  PSD 
Class  II  increments. 

Existing  air  quality  models  were  deemed  inadequate  for  application  to 
the  instantaneous  line  source  emissions  produced  by  intermittent  aircraft 
flights  along  MTRs.  Therefore,  the  MAILS  model  was  created  for  this  unique 
type  of  pollutant  source.  The  model  was  integrated  with  an  aircraft  engine 
emissions  database  to  produce  an  interactive,  user-friendly  modeling  system. 

Validation  of  the  MAILS  model  using  actual  ground-level  concentration 
measurements  was  not  feasible,  because  of  the  lack  of  appropriate  measurements 
and  the  technical  difficulties  in  obtaining  them.  Therefore,  a  performance 
evaluation  of  MAILS  was  accomplished  by  comparing  MAILS  results  with  manually- 
averaged  results  obtained  from  an  existing,  validated,  atmospheric  dispersion 
model  developed  by  the  U.S.  Environmental  Protection  Agency  (EPA).  The  MAILS 
model  test  results  were  in  very  close  agreement  with  the  EPA  model  test 
results. 

The  MAILS  output  provides  estimates  of  air  pollutant  concentrations  from 
proposed  MTRs  and  comparisons  of  the  concentrations  with  existing  PSD  Class  I 
increments.  The  MAILS  model  can  be  run  quickly  and  requires  minimal  user 
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training  to  execute  and  to  interpret  or  apply  results.  The  model  is 
recommended  as  a  planning  and  assessment  tool  for  those  involved  in  the 
environmental  impact  analysis  process  for  proposed  and/or  existing  low- 
altitude  military  airspaces.  An  explanation  of  modeling  techniques  used  to 
calculate  pollutant  concentrations  is  usually  required  in  environmental  impact 
assessments  and  statements.  A  summary  describing  the  MAILS  model,  suitable  for 
attachment  to  such  documents,  is  included  as  an  appendix  to  this  report. 
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This  report  was  prepared  by  the  Oak  Bidge  National  Laboratory,  P.  0.  Box 
2008,  Oak  Ridge,  Tennessee  37831,  under  an  Interagency  Agreement  (lAG)  between 
the  U.S.  Department  of  Energy  (DOE)  and  the  Department  of  the  Air  Force  (DOE 
lAG  Ho.  1A89-1A89-A2 ,  USAF  No.  F88-52),  for  the  Air  Force  Engineering  and 
Services  Center  (AFESC),  Engineering  and  Services  Laboratory,  Tyndall  Air  Force 
Base,  Florida. 
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Officers. 


This  version  of  the  HAILS  software.  Version  3.0,  includes  more  general  menu 
drivers  than  previous  versions.  Version  3.0  and  the  description  herein 
supersedes  all  previous  versions. 


This  report  has  been  reviewed  by  the  Public  Affairs  Office  and  is 
releasable  to  the  National  Technical  Information  Services  (NTIS).  At  NTIS,  it 
will  be  available  to  the  general  public,  including  foreign  nations. 


This  technical  report  has  been  reviewed  and  is  approved  for  publication. 


MICHAEL  T.  MOSS,  Capt,  USAF 
Senior  Research  Meteorologist 


Chief,  Environics  Division 


WAYNE  P.  CHEPREN,  Capt,  USAF 


Chief,  Environmental  Compliance  R&D 


FRANK  P.  GAL^GHER  III,ljfol,  USAF 
Director,  Engineering  and  Services 
Center 


V 

(The  reverse  of  this  page  is  blank.) 


TABLE  OF  CONTENTS 


Section  Page 

I  INTRODUCTION  .  1 

A.  OBJECTIVE .  1 

B.  BACKGROUND  .  1 

C.  SCOPE .  2 

II  MODEL  OVERVIEW  .  5 

A.  DISPERSION  MODEL  .  5 

B.  AIRCRAFT  EMISSIONS  DATABASE  .  5 

C.  MODEL/DATABASE  INTEGRATION  .  6 

D.  COMPUTER  SYSTEM  REQUIREMENTS  .  8 

III  TECHNICAL  DESCRIPTION  .  9 

A.  DISPERSION  MODEL  .  9 

1.  Line  Source  Simulation  .  9 

2.  Concentration  Calculations  .  11 

3.  Averaging  Time  Considerations  .  13 

4.  Model  Input .  14 

5.  Output  Description  .  15 

6.  Model  Performance  Evaluation  .  15 

B.  AIRCRAFT  EMISSIONS  DATABASE  .  21 

1.  Parameters  and  Format .  21 

2.  Emissions  Data  Sources  .  21 

IV  USING  THE  MODEL  AND  DATABASE .  23 

A.  DATABASE  MANIPULATIONS  .  23 

1.  Editing,  Adding,  or  Deleting  Records  .  23 

a.  Emissions  File  .  23 

b.  Emission  Factor  Reference  File  .  26 


TABLE  OF  CONTENTS 
(CONCLUDED) 

Section  Page 

2.  Printing  Database  Sunwnaries .  27 

3.  Reindexing  the  Database .  31 

B.  PERFORMING  A  MODEL  RUN  .  31 

1.  Selecting  Data  for  a  Model  Run  .  31 

2.  Applying  Results  .  35 

C.  EXAMPLE  APPLICATIONS  .  36 

1.  Example  1 .  36 

2.  Example  2 .  41 

V  CONCLUSIONS  .  49 

REFERENCES .  51 

APPENDIX 

A  NAAQS  AND  PSD  CLASS  II  AND  CLASS  I  INCREMENTS .  53 

B  PSD  CLASS  I  AIR  QUALITY  AREAS .  55 

C  MAILS  DISPERSION  MODEL  FORTRAN  CODE  .  65 

D  MAILS  SUMMARY  DESCRIPTION  .  73 


viii 


I 


LIST  OF  FIGURES 

Figure  Page 

1  A  Schematic  of  the  MAILS  Model  and  Database  System  .  .  .  7 

2  Dispersion  Modeling  Concept  .  10 

3  MAILS  Master  Menu  and  Database  Maintenance  Menu  .  24 

4  Single  Emissions  Record  and  Menu  .  25 

5  Emissions  Factor  Document  Reference  Menu  and  Database 

Maintenance  Option  Menu  .  28 

6  Emission  Factor  Reference  Document  Entry  Screen  and 

Text  Input  Screen .  29 

7  Emissions  Data  Summary  Menu  Screen  and  Filter-Setting 

Screen .  30 

8  Emissions  Data  Summary  Screen  after  Filter  Has  Been 

Set  and  Resulting  Summary  Report  .  32 

9  MAILS  Model  Main  Menu  and  Aircraft  Selection  Screen  .  37 

10  Hypothetical  Military  Training  Route  Segments  for 

Example  1 .  38 

11  Aircraft  Data  Entry  Screen .  42 

12  MAILS  SO2  Results  for  Example  1 .  43 

13  Hypothetical  Military  Training  Route  Segments  for 

Example  2 .  44 

14  MAILS  SO2  Results  for  Example  2,  Segment  A  .  47 

15  MAILS  SO2  Results  for  Example  2,  Segment  C  .  48 

B-1  PSD  Class  I  Areas  Designated  Under  the  Clean  Air  Act 

Amendments  of  1977  56 


LIST  OF  TABLES 


Table  Pace 

1  AVERAGING  PERIOD  ADJUSTMENT  FACTORS  .  14 

2  DESCRIPTION  OF  MAILS  MODEL  INPUT  PARAMETERS  .  16 

3  COMPARISON  OF  MAILS  MODEL  AND  ISCST  MODEL  RESULTS  .  20 

4  AIRCRAFT  EMISSIONS  DATABASE  PARAMETERS  .  20 

5  AIRCRAFT  INPUT  DATA  FOR  EXAMPLE  1 .  40 

6  AIRCRAFT  INPUT  DATA  FOR  EXAMPLE  2 .  45 

B-1  PSD  CLASS  I  AREAS  DESIGNATED  UNDER  THE  CLEAN  AIR  ACT 

AMENDMENTS  OF  1977  .  57 


X 


LIST  OF  ABBREVIATIONS  AND  ACRONYMS 


CAA  Clean  Air  Act 

CFR  Code  of  Federal  Pegulations 

CO  carbon  monoxide 

EPA  U.S.  Environmental  Protection  Agency 

GEIS  Generic  Environmental  Impact  Statement 

HC  hydrocarbons 

ISCST 

Industrial  Source  Complex,  Short-Term  (model) 

NAILS 

Multiple-Aircraft  Instantaneous  Line  Source  (model) 

MOA  military  operations  area 

MTR  military  training  route 

NAAQS 

National  Ambient  Air  Quality  Standards 
NO,  nitrogen  dioxide 

PART  particulate  matter 

PC  personal  computer 

PM-IO 

particulate  matter  under  10  microns  in  diameter 

PSD  Prevention  of  Significant  Deterioration  (of  air  quality) 

RA  restricted  area 

SAILS 


Single-Aircraft  Instantaneous  Line  Source  (model) 


sulfur  dioxide 
U.S.  Air  Force 


xi 

(The  Reverse  of  this  page  is  blank) 


SECTION  I 


INTRODUCTION 


A.  OBJECTIVE 

This  report  provides  instructions  for  the  use  of  the  Multiple-Aircraft 
Instantaneous  Line  Source  (MAILS)  atmospheric  dispersion  model.  MAILS  is 
intended  as  a  screening  model  for  prediction  of  air  quality  impacts  of  air 
pollution  emissions  from  low-flying  aircraft  on  military  training  routes 
(MTRs).  An  air  quality  screening  model  is  generally  defined  as  one  that 
provides  worst-case  concentration  estimates  and  can  be  run  quickly  with  a 
minimum  of  user  input.  The  MAILS  model  predictions  of  pollutant 
concentrations  are  considered  to  be  worst-case  estimates  because  of  the 
assumptions  incorporated  in  the  model  design. 

B.  BACKGROUND 

The  impetus  for  development  of  MAILS  was  the  need  to  address  the  air 
quality  impacts  of  low-altitude  MTRs  as  part  of  the  environmental  assessment 
process  required  under  the  National  Environmental  Policy  Act  of  1969  (NEPA, 
Pub.  L.  91-190).  In  accordance  with  NEPA,  the  U.S.  Air  Force  (USAF)  has  begun 
preparation  of  a  Generic  Environmental  Impact  Statement  (GEIS)  addressing  the 
impacts  of  establishing  new  low-altitude  military  airspaces.  Potential  air 
quality  impacts  were  assessed  for  the  preliminary  draft  GEIS  in  part  using  the 
Single-Aircraft  Instantaneous  Line  Source  (SAILS)  dispersion  model  (Reference 
1).  The  SAILS  model  was  later  modified  by  integrating  an  aircraft  emissions 
database  and  adding  the  capability  for  multiple-aircraft  assessment  in  a 
single  model  run.  The  modified  model  was  named  MAILS. 

An  important  finding  from  the  modeling  analysis  conducted  for  the  GEIS  was 
that  the  air  quality  impacts  of  low-flying  military  aircraft  would  be 
negligible  (less  than  5  percent  of  applicable  air  quality  standards)  with 
respect  to  National  Ambient  Air  Quality  Standards  (NAAQS)  and  Prevention  of 
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Significant  Deterioration  (PSD)  Class  II  air  quality  increments,  which  are 
applicable  over  most  areas  of  the  United  States.  The  air  quality  impacts  of 
low-flying  aircraft  were  determined  to  be  potentially  significant  (possibly 
over  5  percent  of  applicable  standards)  only  with  respect  to  the  more 
stringent  limitations  applicable  in  PSD  Class  I  areas  (primarily  national 
parks  and  wilderness  areas  above  certain  sizes).  The  increments  of  pollutant 
concentrations  allowed  in  PSD  Class  I  areas  are  generally  an  order  of 
magnitude  or  more  lower  than  allowed  in  other  areas. 

The  current  NAAQS  and  PSD  Class  II  and  Class  I  increments  are  shown  in 
Appendix  A.  A  listing  and  a  map  of  current  PSD  Class  I  areas  are  shown  in 
Appendix  6.  Standards  occasionally  change  or  are  added,  and  new  areas  can  be 
redesignated  as  PSD  Class  I.  Therefore,  the  air  quality  analyst  is  cautioned 
that  appropriate  state  air  pollution  control  agencies  or  U.S.  Environmental 
Protection  Agency  (EPA)  regional  offices  should  be  consulted  periodically  to 
obtain  any  information  about  new  PSD  Class  I  areas  or  increments. 

An  important  aspect  of  the  PSD  program  is  that  the  allowable  concentration 
increments  apply  only  to  major  stationary  pollutant  sources.  Low-flying 
aircraft,  the  emissions  sources  intended  for  analysis  with  the  MAILS  model, 
are  not  subject  to  any  of  the  regulatory  permitting  and  air  quality  analysis 
requirements  (including  PSD  regulations)  established  pursuant  to  the  Clean  Air 
Act  (CAA)  and  its  amendments.  The  PSD  Class  I  increments  were  established  by 
the  CAA  Amendments  of  1977  (Pub.  L.  95-95)  and  apply  only  to  major  stationary 
sources  of  pollutant  emissions.  However,  to  satisfy  the  air  quality  impact 
analysis  requirements  of  NEPA,  measures  e.stablished  under  the  CAA  (PSD 
increments  and  NAAQS)  were  applied  to  low-flying  aircraft  emissions  to  provide 
a  "yardstick"  for  evaluating  the  significance  of  air  quality  impacts 
(Reference  1). 

C.  SCOPE 

The  types  of  airspaces  considered  in  the  preliminary  draft  GEIS  for  low- 
altitude  flying  operations  included  MTRs,  military  operations  areas  (MOAs), 
and  restricted  areas  (RAs),  among  others.  Existing  atmospheric  dispersion 
models  were  judged  appropriate  for  predicting  ground-level  pollutant 
concentrations  caused  by  low-altitude  flying  operations  within  MOAs  and  RAs. 
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Emissions  from  NOA  and  RA  flying  maneuvers  were  modeled  as  continuous  area 
sources  using  the  ERA  Industrial  Source  Complex,  Short-Term  (ISCST)  dispersion 
model  (Reference  2).  However,  for  analysis  of  MTR  pollutant  impacts,  the 
development  of  a  new  or  modified  model  was  deemed  necessary.  This  decision 
was  based  on  the  characteristics  of  MTR  emissions,  which  consist  of 
intermittent,  essentially  instantaneous  line  sources.  Therefore,  the  MAILS 
model  was  developed  for  application  to  MTRs  or  other  airspaces  or  airspace 
segments  with  consistent  flight  path  orientation  similar  to  MTRs. 
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SECTION  II 


MODEL  OVERVIEW 


A.  DISPERSION  MODEL 

MAILS  Is  an  interactive  air  quality  model  that  can  be  used  for  planning  or 
environmental  assessments  of  MTRs  or  similar  low-altitude  airspaces.  Low- 
altitude  Is  defined  here  as  less  than  3000  feet  above  ground-level.  This 
altitude  Is  consistent  with  the  maximum  altitude  of  operations  addressed 
within  the  context  of  the  GEIS  and  Is  also  roughly  the  upper  limit  of  the 
typical  planetary  boundary  layer,  within  which  atmospheric  turbulence  and 
diffusion  are  largely  dependent  on  the  action  of  wind  and  solar  Insolation  at 
the  earth's  surface.  Most  atmospheric  dispersion  models  are  Intended  for 
application  within  this  layer. 

The  dispersion  algorithms  used  In  MAILS  are  based  on  the  commonly  used 
Gaussian  concentration  distribution  (Reference  3).  The  line  source  emissions 
are  assumed  to  disperse  and  maintain  Gaussian  distributions  in  the  horizontal 
(Y  axis,  perpendicular  to  line  of  release)  and  vertical  directions.  To  obtain 
conservative  concentration  predictions,  the  wind  direction  Is  assumed  to  be 
parallel  to  the  line  of  emissions  release  for  a  period  of  1  hour. 

Concentration  predictions  for  averaging  periods  longer  than  1  hour  are  based 
on  the  frequency  of  expected  flights  as  Input  by  the  model  user  and  on 
empirical  averaging  time  adjustment  factors  to  account  for  meteorological  and 
flight  path  variations  over  these  longer  periods. 

B.  AIRCRAFT  EMISSIONS  DATABASE 

The  aircraft  emissions  database  contains  data  for  roughly  50  military 
aircraft.  Including  some  non-USAF  aircraft.  The  Inclusion  of  non-USAF 
aircraft  in  the  database  Is  necessary  for  assessment  purposes  because  many 
low-altitude  training  routes  are  used  by  aircraft  from  more  than  one  branch  of 
the  military.  The  emissions  database  can  be  easily  updated  through  the  menu- 


driven  MAILS  modeling  system  as  emissions  data  for  new  aircraft,  fuels,  or 
engines  become  available. 

The  database  contains  a  single  record  for  each  aircraft  and  pollutant 
combination.  Since  five  pollutants  are  included  in  the  current  database, 
there  are  approximately  250  total  database  records.  The  pollutants 
represented  in  the  database  are  sulfur  dioxide  (SOj),  nitrogen  dioxide  (NOj), 
particulate  matter  (PART),  carbon  monoxide  (CO),  and  hydrocarbons  (HC).  Only 
the  first  four  pollutants  are  currently  regulated  by  NAAQS,  although  HCs 
contribute  to  photochemical  formation  of  ozone,  for  which  an  NAAQS  has  been 
set.  Only  SOj,  NOj,  and  PART  are  regulated  for  PSD  Class  I  areas.  The  NAAQS 
for  PART  are  now  set  for  particles  under  ten  microns  in  diameter  (PM-10).  The 
PSD  increments  for  PART  are  still  applicable  to  total  particulate  matter, 
although  EPA  has  recently  issued  a  notice  of  proposed  rulemaking  to  replace 
them  with  increments  for  PM-10  (see  Appendix  A). 

C.  MODEL/DATABASE  INTEGRATION 

A  schematic  of  the  NAILS  model  and  database  system  is  shown  in  Figure  1. 
Several  independent  functions  can  be  accomplished  with  this  system:  adding  or 
editing  database  records;  printing  summary  reports  from  the  database;  and 
executing  the  dispersion  model,  with  results  sent  either  to  a  disk  file  or  a 
printer.  The  FoxBASE*  Plus  (Revision  2.0)  database  management  system  was  used 
to  program  the  database  maintenance  features  and  the  model  input  prompts.  The 
executable  (compiled  with  Microsoft*  FORTRAN,  Version  4.1)  MAILS  dispersion 
model  was  incorporated  into  this  structure. 

The  integrated  system  allows  the  user  to  model  the  air  quality  impact  of  a 
single  pollutant  for  several  aircraft  types  in  a  single  model  run.  The  user 
simply  chooses  the  pollutant  to  be  modeled  and  the  number  of  aircraft  types 
for  a  given  route  segment.  The  system  then  displays  a  menu  of  aircraft 
emissions  data  for  selection  and  prompts  for  other  necessary  input,  such  as 
the  frequency  of  flights,  for  the  applicable  averaging  times. 
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MULTIPLE  AIRCRAFT  INSTANTANEOUS  LINE  SOURCE 

(  MAILS  ) 

AIR  QUALITY  MODELING  SYSTEM 


USER  INTERFACE 


Figure  1.  A  Schematic  of  the  MAILS  Model  and  Database  System. 


D.  COMPUTER  SYSTEM  REQUIREMENTS 


The  MAILS  modeling  software  can  be  run  on  systems  meeting  the  following 
requirements: 

1.  an  IBM*-compatible  PC-AT  (or  an  upwardly  compatible  personal  computer 
(PC)  such  as  80386-based  machines)  with  PC-DOS*  or  MS-DOS*  (Version 
3.0  or  later)  operating  system, 

2.  a  minimum  of  360  kB  of  free  memory,  and 

3.  an  Epson*-compatib1e  dot  matrix  printer  (80-column)  or  a  printer  that 
can  emulate  the  Epson*  character  set.  Printer  port  designation  of 
LPTl. 

NOTE;  Although  a  math  coprocessor  (e.g.,  80287,  80387)  is  not  required,  it 
is  highly  recommended  since  it  may  reduce  execution  time  by  a  factor  of  5  or 
6.  If  it  is  known  that  a  math  coprocessor  is  present,  then,  before  executing 
MAILS,  at  the  DOS  prompt,  the  user  should  type  COPY  SAILW.EXE  SAILS.EXE.  If 
tt.e  math  coprocessor  is  not  present,  the  user  should,  instead,  type  COPY 
SAILWO.EXE  SAILS.EXE.  If  it  is  not  known  whether  the  coprocessor  is  present, 
assume  it  is  not. 

The  MAILS  software  should  be  copied  to  and  run  from  a  hard  disk.  It  is 
advisable  to  create  a  separate  subdirectory  and  copy  all  files  on  the  MAILS 
diskettes  to  this  subdirectory  (see  your  operating  system  user's  documentation 
or  print  the  README.DOC  file  from  ti »  MAILS  diskettes  for  instructions). 
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SECTION  III 


TECHNICAL  DESCRIPTION 


A.  DISPERSION  MODEL 

This  section  describes  the  design  of  the  MAILS  atmospheric  dispersion 
model,  the  formats  of  model  input  and  output,  and  a  model  performance 
evaluation.  A  listing  of  the  MAILS  dispersion  model  FORTRAN  computer  code  is 
provided  in  Appendix  C. 

1.  Line  Source  Simulation 

The  MAILS  dispersion  model  calculates  worst-case  (maximum  potential) 
1-hour  concentrations  for  linear  flight  paths  assumed  to  be  parallel  to  the 
wind  direction.  For  a  crosswind  situation,  a  ground-level  receptor  would  be 
impacted  for  only  a  brief  time  after  each  aircraft  pass,  resulting  in  much 
lower  1-hour  concentrations  than  for  the  parallel  wind  case.  Concentrations 
for  averaging  periods  longer  than  1  hour  are  obtained  by  accounting  for  (1) 
the  numbers  of  flights  over  these  longer  periods  and  (2)  variations  in 
meteorology  and  flight  paths. 

Emissions  from  a  single  flight  path  can  be  characterized  as  an 
essentially  instantaneous,  infinite  line  source.  The  MAILS  model  divides  the 
instantaneous  line  source  into  puffs,  each  having  a  pollutant  mass  equal  to 
the  emissions  contained  in  a  100-meter-long  segment  of  flight  path.  This 
dispersion  modeling  concept  is  illustrated  in  Figure  2.  For  calculating  the 
impacts  of  multiple  flights  along  the  same  route  path  over  periods  of  1  hour 
or  longer,  the  model  simply  sums  the  calculated  1-hour  contributions  from 
individual  flights  and  applies  averaging  period  adjustment  factors  as 
described  later  in  this  section. 

For  a  single  aircraft  pass,  the  MAILS  model  calculates  the 
concentrations  for  an  array  of  wind  speed/atmospheric  stability  conditions  and 
then  selects  the  maximum  concentration  for  any  of  these  meteorological 
conditions.  These  49  conditions  are  identical  to  those  used  in  the  PTPLU 


9 


TIME  =  ti  WIND  DIRECTION 


Figure  2.  Dispersion  Modeling  Concept. 


model  (Reference  4),  which  Is  an  EPA-approved  screening  model  for  single 
continuous  point  sources.  The  MAILS  dispersion  model  assumes  wind  speed  and 
direction  to  be  constant  In  space  and  time. 

Concentrations  are  computed  for  each  wind  speed/stability  condition  at 
a  ground-level,  plume  centerline  location  (receptor)  by  summing  the 
contributions  of  all  the  puffs  that  pass  the  location  In  a  1-hour  period.  The 
contributions  to  total  exposure  from  puffs  more  than  1  hour  of  travel  time 
away  are  expected  to  be  small  for  most  meteorological  conditions  and  were  not 
considered  In  obtaining  maximum  concentration  estimates  for  periods  longer 
than  1  hour.  Furthermore,  because  low-altitude  wind  directions  are  unlikely 
to  remain  constant  for  an  hour  or  more,  the  same  receptor  would  probably  not 
be  exposed  to  Impacts  from  the  centerline  of  the  puffs/plume  for  longer 
periods.  Thus,  the  concentrations  calculated  by  the  model  for  a  single  linear 
flight  path  are  considered  to  be  conservative  (upper-bound)  1-hour  estimates. 

2.  Concentration  Calculations 

The  MAILS  model  calculates  concentrations  by  summing  the  exposures 
(product  of  concentration  and  time)  from  Individual  puffs  that  pass  the 
receptor  In  a  1-hour  period  and  by  converting  the  total  exposure  to  a  1-hour 
concentration.  The  mass  of  each  puff  Is  assumed  to  be  distributed  according 
to  a  Gaussian  shape  In  the  vertical  (z)  and  horizontal  (x  and  y)  directions. 
The  exposure  V  (9  sec)  from  an  Individual  puff  is  given  by  the  following 
equation  taken  from  Slade  (Reference  5,  p.  115): 


where  Q  Is  the  puff  mass  (grams),  Is  the  standard  deviation  (m)  of  the 
horizontal  puff  concentration  distribution,  Is  the  standard  deviation  (m) 
of  the  vertical  puff  concentration  distribution,  u  Is  the  average  wind  speed 
(m/s),  h  Is  the  emissions  release  height  (m)  above  ground-level,  and  y  is  the 
perpendicular  distance  (m)  from  the  receptor  to  the  plume  centerline.  Because 
only  the  concentration  on  the  plume  centerline  (at  ground  level)  is  calculated 
by  the  MAILS  model,  y  Is  zero  and  the  first  term  in  the  exponent  drops  out. 
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The  preceding  equation  is  appropriate  for  a  situation  with  essentially 
unlimited  vertical  atmospheric  mixing.  Where  limited  mixing  is  important,  an 
expansion  of  the  exponential  term  is  necessary  to  account  for  multiple 
reflections  of  the  Gaussian  vertical  puff  profile  off  the  mixing  height 
inversion  and  ground  surfaces.  In  such  cases  the  exponential  term  is  expanded 
in  the  same  manner  as  done  by  the  ISCST  model  (Reference  2).  The  ground- 
level  plume-centerline  exposure  is  then  given  by 


where  H  is  the  surface  mixing  or  inversion  height. 

At  downwind  distances  where  the  ratio  of  is  greater  than  1.6, 

the  puff  is  essentially  uniformly  mixed  in  the  vertical  direction.  The  plume 
centerline  exposure  from  a  single  puff  is  then  given  by  the  equation 


Q 

^  =  - z — 

ir(TyU  H 


(3) 


The  dispersion  coefficients  ((Tyand  a^)  used  in  the  MAILS  model  are 
identical  to  those  used  by  the  EPA  ISCST  dispersion  model  (Reference  2}  when 
the  ISCST  model  is  applied  to  point  source  plumes  without  downwash- induced 
dispersion  enhancement,  in  the  "rural"  mode.  The  use  of  the  rural  dispersion 
coefficients  is  considered  appropriate  for  this  application  because  the  MTRs 
to  which  the  MAILS  model  is  applied  would  not  cross  over  urban  areas.  It 
might  be  argued  that  dispersion  coefficients  for  puffs,  which  are  somewhat 
smaller  than  for  continuous  plumes,  should  be  used  in  the  MAILS  model. 
However,  it  is  appropriate  to  use  dispersion  coefficients  for  plumes  rather 
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than  for  puffs  because  the  puffs  are  used  by  the  model  to  simulate  a 
continuous  plume  produced  by  the  essentially  instantaneous  line  source 
emissions. 

3.  Averaging  Period  Adjustments 

Air  quality  standards  are  typically  established  for  averaging  periods 
ranging  from  1  hour  to  1  year  (see  Appendix  A).  The  HAILS  model  estimates 
maximum  concentrations  for  various  averaging  periods  by  summing  the 
contributions  of  individual  flights  that  would  occur  over  the  period  of 
Interest  and  by  applying  empirically  derived  averaging  period  adjustment 
factors.  The  maximum  concentration  for  a  particular  averaging  period  is  given 
by 

n 

i=l 

where  X  is  the  maximum  concentration  (m1crograms/mi^),  n  Is  the  number  of 
different  aircraft/altitude  combinations  being  modeled,  C,  Is  the  maximum  1- 
hour  concentration  contribution  (mlcrograms/m’)  from  a  single  overflight  for 
one  aircraft/altitude  combination,  F,  Is  the  frequency  of  overflights  for  one 
aircraft/altitude  combination  during  the  period  of  Interest,  P  is  the  period 
of  interest  (hours),  and  A  is  the  averaging  period  adjustment  factor  for  the 
period  P.  The  averaging  period  adjustment  factors  account  for  the  variations 
in  meteorological  conditions  and  in  the  horizontal  and  vertical  position  of 
subsequent  flights  with  respect  to  the  route  centerline.  These  factors  are 
subjective  estimates  based  on  dispersion  modeling  experience  and  based  on  EPA 
recommendations  (Reference  6)  for  other  types  of  pollutant  sources.  The 
averaging  period  adjustment  factors  used  for  various  periods  are  shown  in 
Table  1. 

The  3-hour,  8-hour  and  24-hour  factors  are  somewhat  lower  than  similar 
factors  recommended  by  EPA  for  relating  maximum  1-hour  concentrations  to 
longer  period  maximum  concentrations  associated  with  longer  periods  for 
continuous  elevated  point  sources  (Reference  6).  The  EPA  factors  are  intended 
to  account  for  variations  in  meteorological  conditions  that  cause  longer-term 
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TABLE  1.  AVERAGING  PERIOD  ADJUSTNENT  FACTORS 


1  hour 

1.00 

3  hour 

0.50 

8  hour 

0.33 

24  hour 

0.25 

Annual 

0.10 

concentrations  to  be  lower  than  1-hour  concentrations.  However,  the  EPA 
factors  were  developed  for  continuous  stationary  sources,  while  the  factors  in 
Table  1  were  developed  for  an  intermittent  source,  for  which  each  successive 
emission  varies  in  both  horizontal  and  vertical  position  relative  to  a  fixed 
ground-level  receptor.  The  variation  in  horizontal  and  vertical  position  of 
the  successive  line  sources  creates  greater  effective  dispersion  of  the 
emissions  as  averaging  time  increases,  thereby  lowering  further  the 
appropriate  averaging  period  adjustment  factors. 

The  Table  1  averaging  period  adjustment  factors  are  not  appropriate  in 
situations  involving  extremely  low  flight  frequencies.  For  example,  if  a 
maximum  of  one  flight  is  expected  in  a  3-hour  period,  flight  path  variation 
does  not  apply,  and  variations  in  meteorological  conditions  over  periods 
greater  than  1  hour  are  not  relevant.  In  such  a  case  it  would  be  appropriate 
to  set  A  equal  to  1.0.  However,  because  routes  with  very  low  flight 
frequencies  would  generate  negligible  impacts  with  respect  to  PSD  Class  I 
increments,  this  model  limitation  is  not  a  problem  for  NEPA  analyses  or  other 
applications  where  PSD  Class  I  increments  are  used  as  the  measure  of  impact 
intensity. 

4.  Model  Input 


This  section  describes  the  input  parameters  required  by  the  MAILS 
dispersion  model.  Some  of  these  parameters  can  be  selected  from  the  emissions 


database  or,  at  the  user's  option,  edited  upon  selection  of  an  aircraft 
emissions  record  for  a  problem  run  without  actual  modification  of  the  database 
file.  Other  parameters  must  be  input  by  the  user,  based  on  the 
characteristics  of  the  route  segment  to  be  modeled. 

Table  2  lists  the  model  input  parameters  in  the  sequence  in  which  the 
user  is  prompted  by  the  model.  Where  the  parameter  type  is  shown  as  "0,"  a 
default  parameter  is  retrieved  from  the  database  upon  user  selection  of  an 
aircraft-pollutant  record  from  a  menu  and  can  be  modified  before  input  to  the 
model  at  the  option  of  the  user.  Where  the  parameter  type  is  "U,"  the  value 
for  the  parameter  is  a  user  input. 

5.  Output  Description 

The  MAILS  model  output  is  sent  either  to  a  file  named  by  the  user 
during  model  input,  or  directly  to  the  printer.  If  the  printer  option  is 
chosen,  the  default  output  device  designation  is  LPTl.  As  a  backup,  in  case 
the  printer  cannot  be  accessed,  the  output  is  also  sent  to  a  file  named 
MODEL. PRT  when  the  printer  option  is  chosen.  Some  examples  of  MAILS  model 
output  are  shown  later  in  Figures  12,  14,  and  15.  The  model  output  consists 
of  summary  tables  for  each  of  the  applicable  averaging  periods.  Each  table 
shows  all  model  input;  the  calculated  concentration  contributions  by  aircraft 
type/altitude  combination;  total  concentration  for  the  MTR  segment  or  set  of 
concurrent  segments;  and,  if  the  standard  run  mode  was  selected,  the 
calculated  percentage  of  NAAQS  and  PSD  Class  I  increment  consumed  by  the 
modeled  MTR  segment(s).  If  the  calculated  percentage  of  Class  I  increment 
consumption  is  less  than  5  percent  for  a  given  pollutant,  the  MTR  segment(s) 
is  considered  to  have  a  negligible  impact  with  respect  to  that  pollutant  in 
the  Class  I  area.  A  more  detailed  discussion  of  the  application  of  model 
output  is  contained  in  Sect.  IV  (B). 

6.  Model  Performance  Evaluation 

For  many  types  of  dispersion  models,  at  least  some  field  measurements 
are  available  for  an  objective  assessment  of  model  performance.  However,  for 
some  types  of  models  the  necessary  field  data  do  not  exist  and  would  be 
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TABLE  2.  DESCRIPTION  OF  HAILS  MODEL  INPUT  PARAMETERS 


Parameter _ Type 


Title  U 

Run  mode  U 


Pollutant  U 


Output  U 

option 


Number  of  U 

aircraft 


Aircraft  D 


_ Description _ 

User's  description  of  the  model  run  (e.g., 
airspace  description,  date  of  run). 

User  can  select  standard  (S)  or  nonstandard 
(N)  mode.  Choice  of  the  standard  mode 
results  in  prompts  for  flight  frequency  (see 
“Frequency"  parameter  below)  only  for  the 
averaging  periods  and  pollutants  with 
applicable  air  quality  standards  (see 
Appendix  A).  Choice  of  the  nonstandard  mode 
will  result  in  prompts  for  frequency  for  i- 
hour,  3-hour,  8-hour,  24-hour,  and  annual 
periods  and  allows  modeling  of  other 
pollutants. 

The  database  contains  emissions  data  for  six 
pollutants.  A  single  pollutant  (CO,  HC,  NOg, 
PART,  or  SOg)  is  selected  for  a  model  run. 
Only  three  pollutants  (SO,,  NO,,  and  PART) 
must  be  considered  for  a  PSD  Class  I 
analysis.  The  user  is  allowed  to  choose 
"OTHER"  only  if  "N"  was  selected  for  Run  mode 
(above) . 

The  user  enters  P  to  send  output  directly  to 
the  printer  (a  backup  file  named  MODEL. PRT  is 
also  written  in  this  case)  or  F  to  send 
output  to  a  file.  If  F  is  selected,  the  user 
is  prompted  for  the  file  name  and  that  file 
name  is  then  suffixed  with  .DAT. 

This  value  is  actually  the  number  of 
aircraft/altitude  combinations  for  the 
military  training  route  (MTR)  segment  to  be 
modeled.  In  nearly  all  situations,  a  given 
aircraft  type  is  assigned  only  one  minimum 
altitude  for  a  given  route  or  set  of 
concurrent  route  segments. 

The  aircraft  designation,  for  example,  "B- 
52G,"  in  database  is  selected  by  user. 
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TABLE  2.  DESCRIPTION  OF  NAILS  NODEL  INPUT  PARANETERS 

(CONTINUED) 


Parameter _ Type 

Airspeed  D 

Altitude  U 

Nixing  U 

height 

Emission  D 

rate 

Frequency  U 


_ Description _ 

Average  airspeed  (miles/hour)  of  the 
particular  aircraft  on  MTR  operations. 

The  minimum  altitude  (feet)  allowed  for  the 
particular  aircraft  on  the  MTR  segment.  If 
the  minimum  altitude  is  less  than  200  feet,  a 
value  of  200  feet  is  recommended  for  input 
unless  the  aircraft  consistently  flies  below 
this  altitude. 

The  height  (in  feet)  of  the  surface-based 
temperature  inversion.  Use  of  the  default 
value  of  5000  feet  (essentially  unlimited 
vertical  mixing)  is  recommended  for  most 
applications. 

Total  emission  rate  (Ib/hour)  of  a  single 
aircraft  (all  engines)  for  the  selected 
pollutant. 

Number  of  passes  by  the  designated  aircraft 
type  on  the  route  segment  for  each  applicable 
averaging  period,  as  determined  by  run  mode 
and  pollutant  selection.  For  example,  if  the 
user  selects  "standard"  for  run  mode  and 
"NOj"  for  pollutant,  the  user  will  be 
prompted  only  for  an  annual  frequency,  since 
only  annual  standards  exist  for  NOj.  The 
user  is  cautioned  that  this  is  not  always  the 
same  as  the  number  of  sorties  (flights  or 
missions).  For  example,  some  MTRs  have 
"racetrack"  portions  where  one  aircraft  may 
make  several  passes  during  a  given  sortie. 
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prohibitively  expensive,  if  not  impossible,  to  gather.  This  is  the  '  .se  with 
the  MAILS  dispersion  model. 

Some  limited  field  experiments  have  been  conducted  to  measure  ground- 
level  concentrations  produced  by  elevated  instantaneous  line  sources.  Such 
sources  were  simulated  by  tracers  released  from  low-flying  aircraft  (Reference 
5,  pp.  170-173).  However,  these  experiments  were  configured  so  that  the  line 
of  tracer  release  was  perpendicular  to  the  wind  direction.  Vertically  spaced 
samplers  mounted  on  a  tower  and  a  line  of  ground-based  samplers  downwind  of 
the  tracer  releases  recorded  the  concentrations  produced  by  the  tracer  as  it 
dispersed.  The  MAILS  model  dispersion  algorithm  simulates  worst-case,  short¬ 
term  concentrations  produced  by  an  essentially  instantaneous  elevated  release 
that  is  parallel  to  the  wind  direction.  Thus,  as  shown  in  Figure  2,  a  fixed 
ground-level  receptor  under  the  plume  centerline  would  be  exposed  to  emissions 
for  a  much  longer  period,  as  "puffs"  originating  from  successively  farther 
upwind  locations  continue  to  pass  over  the  receptor.  The  difficulty  in 
conducting  a  field  experiment  to  measure  such  impacts  would  be  in  creating  a 
tracer  release  so  that  the  centerline  of  the  emissions  would  be  exactly 
parallel  to  the  wind  direction  for  several  miles  upwind  of,  and  in  line  with, 
a  ground-based  monitor  (receptor). 

Because  of  the  lack  of  appropriate  field  data  for  model  comparison,  an 
alternative  model  evaluation  approach  was  utilized.  This  approach,  in  accord 
with  the  general  guidance  provided  by  EPA  (Reference  7),  consisted  of 
comparing  the  results  of  the  MAILS  model  with  the  results  of  the  IS^ST  model 
(Reference  2),  which  has  already  been  the  subject  of  a  performance  evaluation 
and  which  can  be  configured  to  simulate  the  type  of  source  to  which  the  MAILS 
model  will  be  applied. 

The  following  inputs  and  procedures  were  used  to  obtain  ISCST  results 
that  could  be  compared  to  those  produced  by  the  MAILS  model. 

a.  Both  models  were  run  for  a  single  source. 

b.  ISCST  comparison  results  were  obtained  using  the  "rural" 
dispersion  mode  and  a  wind  speed  of  10  m/s,  under  two 
atmospheric  stability  categories  [Pasquill  C  and  D  (see 
Reference  3)],  under  both  limited  and  unlimited  mixing 


situations.  Thus,  four  meteorological  cases  were  generated  for 
comparison. 

c.  An  ISCST  emission  rate  was  calculated,  based  on  the  10  m/s  wind 
speed,  to  provide  the  same  linear  emissions  density  produced  by 
MAILS  Input  for  the  comparison  run. 

d.  The  ISCST  structural  downwash  options  were  not  used,  so  that  the 
Gaussian  plume  dispersion  parameters  would  be  the  same  as  used 
by  the  MAILS  model. 

e.  A  300-foot  (91.5-m)  stack  source,  with  no  plume  rise  (near-zero 
exit  velocity  and  stack  diameter)  was  simulated  with  ISCST  to 
produce  a  plume  of  constant  height. 

f.  Receptors  were  placed  every  500  meters  In  the  downwind  direction 
for  the  ISCST  run,  out  to  a  distance  of  36  km,  which  would  be 
the  distance  of  puff  transport  In  1  hour  with  the  10  m/s  wind 
speed. 

g.  The  72  concentration  results  produced  by  ISCST  at  the  receptors 
described  In  Item  f  were  added  and  divided  by  72  to  yield  values 
for  comparison  with  the  MAILS  output. 

Use  of  the  ISCST  results  In  the  preceding  manner  produces  an  average 
concentration  at  a  receptor  that  is  moving  along  the  centerline  of  a 
dispersing  plume  of  constant  linear  density  for  a  period  of  1  hour.  The  MAILS 
model,  on  the  other  hand,  simulates  the  concentration  at  a  fixed  receptor 
under  the  centerline  of  a  dispersing  plume  of  constant  linear  density,  which 
Is  assumed  to  move  over  the  receptor  for  1  hour.  Thus,  these  converse 
simulations  should  produce  comparable  results,  and  the  verification  Is  shown 
in  Table  3.  The  differences  between  the  two  model  results  are  less  than  5 
percent  for  all  cases.  Indicating  that  the  MAILS  model  performs  comparably 
with  an  existing  validated  dispersion  model.  This  result  Is  not  unexpected, 
because  Identical  dispersion  coefficients  are  used  In  the  MAILS  model  and  In 
the  ISCST  model  as  executed  for  this  comparison.  Use  of  the  ISCST  model  for 
routine  air  quality  analyses  of  MTRs  Is  undesirable  for  several  reasons;  the 
ISCST  model  Is  non-InteractIve,  ISCST  requires  more  detailed  and  complex  input 
than  does  MAILS,  ISCST  Is  a  much  larger  and  therefore  less  portable  program 
code  as  compared  with  MAILS,  and  manual  averaging  of  ISCST  concentration 


19 


TABLE  3.  COMPARISON  OF  MULTIPLE  AIRCRAFT  INSTANTANEOUS  LINE  SOURCE 

(NAILS)  AND  INDUSTRIAL  SOURCE  COMPLEX, 

SNORT-TERN  (ISCST)  MODEL  RESULTS 


Stability 

Mixing  height 
(m) 

MAILS 

concentration 

(micrograms/tn^) 

ISCST 

concentration 

(micrograms/m^) 

C 

122 

5.072 

5.069 

C 

5000 

1.781 

1.851 

D 

122 

5.883 

5.899 

D 

5000 

3.493 

3.501 

TABLE  4. 

AIRCRAFT  EMISSIONS  DATABASE  PARAMETERS 

Parameter 

Type 

Field 

width 

Decimal 

places 

Units 

(if  applicable) 

Aircraft  type 

Character 

8 

Average  speed 

Number 

3 

0 

mph 

Number  of  engines 

Number 

1 

0 

Fuel  rate  per  engine 

Number 

8 

2 

1000  lb  fuel/h 

Emission  factor 

Number 

8 

2 

lb/1000  lb  fuel 

Emiss.  factor  flag 

Character 

1 

Pollutant 

Character 

4 

Emission  rate 

Number 

8 

2 

Ib/h 

Emission  density 

Number 

6 

4 

Ib/mile 
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estimates  at  receptors  was  necessary  to  produce  the  values  for 
comparison. 

B.  AIRCRAFT  EMISSIONS  DATABASE 

1.  Parameters  and  Format 

The  aircraft  emissions  database  contains  one  record  for  each 
aircraft/pollutant  combination.  The  parameters,  with  their  formats  and  units, 
are  listed  in  Table  4.  Most  of  the  parameters  in  Table  4  are  self- 
explanatory.  The  emission  rate  is  a  product  of  the  number  of  engines,  the 
fuel  rate  per  engine,  and  the  emission  factor,  which  is  given  in  pounds  of 
pollutant  per  thousand  pounds  of  fuel.  The  emission  factor  flag  is  a  single 
character  that  is  referenced  to  a  separate  file  containing  text  to  explain  the 
source  of  each  emission  factor  in  the  database. 

The  emissions  density  is  not  used  by  the  model  directly  but  was 
included  in  the  database  to  allow  the  user  to  screen  the  emissions  data 
visually  to  determine  potential  worst-case  MTR  segments.  If  the  operating 
altitudes  of  two  different  aircraft  are  equal,  the  emissions  densities  for 
these  aircraft  will  provide  a  comparison  of  relative  ground-level 
concentration  caused  by  the  aircraft  (i.e.,  if  Aircraft  A  has  an  emissions 
density  twice  that  of  Aircraft  B,  the  ground-level  concentration  caused  by 
Aircraft  A  will  be  twice  that  of  Aircraft  B).  When  a  record  is  added  to  the 
database,  the  emissions  density  is  calculated  automatically  by  the  database 
software  by  dividing  the  emission  rate  by  the  average  airspeed. 

2.  Emissions  Data  Sources 

Most  of  the  aircraft  emissions  data  are  based  on  emission  factors  and 
fuel  rates  drawn  from  Seitchek  (Reference  8).  These  data  have  been 
supplemented  with  data  from  other  sources  (References  9-11).  The  SO2  emission 
rates  for  all  aircraft  engines  were  calculated  based  on  a  fuel  sulfur  content 
of  0.05  percent,  which  is  on  the  high  end  of  the  numerous  aviation  fuel  test 
results  reported  by  Shelton  (Reference  9).  A  separate  emission  factor 
reference  file  contains  a  brief  description  of  the  data  source  or  the 
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assumptions  used  in  obtaining  or  calculating  each  emission  factor.  Each  entry 
In  the  emission  factor  reference  file  is  Identified  by  a  one  character  "flag" 
that  corresponds  with  the  emission  factor  flag  used  in  the  main  database. 

Because  air  pollutant  emission  rates  vary  substantially  with  engine 
throttle  setting,  it  is  important  to  choose  data  representative  of  the 
throttle  settings  used  for  typical  low-altitude  flight  operations.  The 
emission  rates  in  the  database  are  representative  of  the  "intermediate"  engine 
operating  mode  as  used  in  Reference  8.  The  intermediate  mode  was  considered 
to  be  the  most  applicable  to  the  type  of  cruising  operations  that  would  occur 
on  most  MTRs.  However,  data  from  other  sources  were  not  generally  categorized 
using  the  same  terminology  used  in  Reference  8.  For  these  other  sources,  a 
data  point  was  chosen  that  most  closely  corresponded  to  approximately  75 
percent  of  the  maximum  throttle  or  fuel  flow  rate  under  "normal"  engine 
operating  conditions,  apart  from  afterburner  or  other  special  modes  not 
intended  for  continuous  operation. 

The  emissions  data  from  Seitchek  and  other  sources  are  used  in  a 
conservative  manner  in  the  MAILS  model.  SOg  emissions  are  based  on  complete 
conversion  of  the  fuel  sulfur  to  SOj.  NOg  emissions  are  based  on  the 
assumption  that  all  nitrogen  oxides  (NO  and  NOj)  are  in  the  form  of  NOg. 
Finally,  total  particulate  matter  emissions  are  used  to  produce  results  for 
comparison  against  standards  for  particulate  matter  under  10  microns  in 
diameter  (PM-10).  These  conservative  assumptions  would  typically  result  in 
SO2  concentrations  that  are  several  percent  higher  than  actual  concentrations, 
and  NOj  and  PM-10  concentrations  up  to  roughly  a  factor  of  two  higher  than 
actual  concentrations. 


22 


SECTION  IV 


USING  THE  MODEL  AND  DATABASE 
A.  DATABASE  MANIPULATIONS 

Most  of  the  aspects  of  database  manipulation  are  easily  understood  by 
moving  through  the  menu  screens  and,  if  necessary,  using  the  "Help"  options 
along  the  way.  However,  explanations  of  some  of  the  common  types  of 
manipulations  are  given  in  the  following  sections. 

1.  Editing,  Adding,  or  Deleting  Records 
a.  Emissions  File 

Emissions  file  records  can  be  easily  modified,  added,  or  deleted 
by  selecting  the  "Database  Maintenance"  feature  from  the  MAILS  master  menu  and 
the  "Revise  Emissions  Database"  feature  from  the  DATABASE  MAINTENANCE  menu 
that  follows.  These  menus  are  shown  in  Figure  3.  The  program  then  prompts 
the  user  with  a  screen  listing  the  various  database  parameters  for  the  first 
record  in  the  database.  An  example  of  this  screen  is  shown  in  Figure  4.  The 
fields  on  this  first  screen  cannot  be  modified  at  this  point;  they  are 
intended  only  to  give  the  user  an  indication  of  the  type  of  data  in  the 
database.  The  user  can  use  the  menu  features  at  the  bottom  of  the  screen  to 
scan  the  existing  data  or  to  add,  copy,  modify,  or  delete  records.  For 
example,  selecting  the  "add"  feature  clears  the  aircraft  identification  field 
and  most  of  the  other  data  fields,  except  for  some  default  values.  The  user 
can  then  enter  the  desired  data  in  the  empty  fields,  concluding  with  the 
emission  factor  flag.  The  emission  rate  is  not  entered  by  the  user  when 
adding  a  record;  the  database  computes  this  value  automatically  by  taking  the 
product  of  the  fuel  rate  per  engine  (1000  Ib/h),  the  number  of  engines,  and 
the  emission  factor  (pounds  of  pollutant/1000  lb  fuel). 
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Air  Quality  Modeling  System 


Run  MAILS  Model 
DataBase  Maintenance 
Introductory  Screen 
Quit  ... 


Air  Quality  Modeling  System 


DATABASE  MAINTENANCE 


Revise  Emissions  DataBase 
Emissions  Flag  Ref. 

Manual  Reindex  Function 
Print  Emissions  Database  Summary 
Help  Screen  Review 
Quit  ... 


*  Maint  Menu  * 


Figure  3.  HAILS  Master  Menu  (top)  and  Database  Maintenance  Menu  (bottom). 
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Air  Quality  Modeling  System 


Aircraft  Emissions  Data 


AIRCRAFT  IDENTIFICATION 

. ---> 

AlO 

AIRSPEED  . 

. > 

405  (  MPH  ) 

FUEL  RATE  PER  ENGINE  - 

. > 

0.92  (1000  Lb/Hr) 

EMISSION  FACTOR  . 

. > 

16.30  (Lb/lOOOLb) 

TYPE  or  POLLUTANT  . 

. > 

CO 

NUMBER  OF  ENGINES  . 

->  ;  2 

EMISSION  FLAG  --->  :  B 

EMISSION  RATE  . 

. > 

29.99  (Lb/Hr) 

Ret/Beg/End/Nxt/Prv/Skp/Mod/Add/Cpy/Del/Lst/Fi 1 t/Tal ly/Hel p/Vi ew/Qui t 
Retrieve  a  record  by  its  key  fields 


Figure  4.  Single  Emissions  Record  and  Nenu. 
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Another  way  to  add  a  record  is  to  use  the  "copy"  feature.  This 
is  especially  useful  when  adding  a  new  record  to  the  database  that  is  very 
similar  to  an  existing  record.  In  this  case  the  user  scrolls  through  the 
database  to  pull  up  the  similar  record  on  the  screen.  Selecting  "Copy"  will 
cause  these  data  to  be  copied  to  an  identical  data  entry  screen  and  allow  the 
user  to  make  changes  in  any  of  the  parameters.  When  the  cursor  is  moved  to  a 
point  before  the  first  data  field  or  past  the  last  data  field,  the  software 
asks  the  user  if  the  record  is  to  be  saved;  if  the  user  answers  "Y,"  the 
record  is  added  to  the  permanent  database. 

An  easy  way  to  scan  the  database  is  to  select  the  "Lst"  feature 
and  then  use  the  up/down  arrow  keys  or  the  page  up/page  down  keys  to  browse 
the  database  quickly.  Once  a  record  is  highlighted  with  this  feature 
pressing  the  "Esc"  key  will  cause  this  record  to  be  retrieved  to  the  single 
record  display  screen  for  possible  manipulation. 

There  are  a  number  of  other  useful  features  in  the  menu  at  the 
bottom  of  the  single  record  display  screen  (Figure  4).  A  one-line  explanation 
of  each  function  is  given  as  the  function  is  highlighted  by  using  the  arrow 
keys.  Also,  if  the  user  wants  to  select  a  particular  function  from  this  menu, 
this  can  be  done  either  by  highlighting  the  function  and  pressing  "Enter"  or 
by  simply  typing  the  first  letter  of  the  desired  function. 

The  "View"  function  allows  the  user  to  see  the  reference 
documentation  for  the  emission  factor.  Documentation  of  the  source  of  each 
emission  factor  in  the  database  is  contained  in  the  emission  factor  reference 
file.  Editing  and  printing  of  the  emission  factor  reference  file  is  described 
in  the  following  section. 

b.  Emission  Factor  Reference  File 

This  file  documents  the  sources  of  various  emission  factor  data. 
If  the  emissions  data  file  is  to  be  revised,  the  user  should  first  print  a 
summary  of  the  emission  factor  reference  file  and  determine  if  one  of  the 
existing  references  in  this  file  is  applicable  as  a  source  for  a  new  or 
revised  emission  file  data.  If  one  of  the  existing  references  is  applicable, 
the  user  should  enter  the  corresponding  single-character  flag  when  the 
emissions  record  is  added  or  revised.  If  none  of  the  existing  emission  factor 
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references  are  appropriate  for  the  new  data,  the  user  should  add  a  new  flag 
and  explanation  to  the  emission  factor  reference  file.  This  can  be  done  by 
selecting  the  "Emission  Flag  Ref."  feature  from  the  DATABASE  MAINTENANCE 
screen  (Figure  3).  The  user  will  then  be  prompted  by  the  Emissions  Factor 
Document  Reference  menu  (Figure  5,  top),  from  which  the  user  can  print  a 
summary  of  the  reference  file  (Summary  Report  option)  or  can  revise  the 
database  (Database  Maintenance  option).  If  the  latter  is  chosen,  the  screen 
shown  at  the  bottom  of  Figure  5  is  shown. 

Adding,  modifying,  or  deleting  new  flags  is  done  in  a  similar 
manner  as  for  the  emissions  data  file  (see  Sect.  IV  A.l.a.).  Adding  or 
editing  the  text  that  accompanies  the  flag  is  done  differently.  For  example, 
to  add  a  flag  and  accompanying  reference  text,  the  user  selects  the  "Add" 
feature  (see  bottom  of  Figure  5)  and  the  software  blanks  out  the  Emissions 
Reference  Flag  field.  After  the  desired  flag  is  entered,  the  user  is  asked  if 
he  would  like  to  edit  the  document  area.  If  the  answer  is  yes,  the  screen  at 
the  top  of  Figure  6  is  shown.  The  user  must  then  press  the  "Ctrl"  and  "Home" 
keys  to  display  the  text  entry  screen,  shown  at  the  bottom  of  Figure  6.  After 
the  desired  text  is  entered,  the  user  must  press  the  "Ctrl"  and  "W"  keys  to 
write  the  text  to  the  database  record  associated  with  the  selected  flag. 

If  the  "Help"  feature  is  chosen  from  the  menu  shown  at  the 
bottom  of  Figure  5,  a  brief  summary  of  the  options  and  keystrokes  needed  to 
edit  the  Emission  Factor  Reference  Database  is  provided. 

2.  Printing  Database  Summaries 

To  obtain  a  hard  copy  listing  of  emissions  data,  the  user  must  select 
the  "Print  Emissions  Database  Summary"  option  (Figure  3)  from  the  DATABASE 
MAINTENANCE  menu.  The  user  can  then  browse  the  data  using  the  "Lst"  and  other 
options  that  will  appear  as  shown  at  the  top  screen  in  Figure  7.  If  "Go"  is 

selected  at  this  point,  a  tabular  summary  will  be  produced  of  all  records  in 

the  emissions  database.  However,  the  user  can  print  a  summary  of  a  selected 

set  of  records  by  first  selecting  the  "Filt"  or  filter  option.  A  screen  will 

be  displayed  with  certain  database  fields  empty,  but  highlighted.  The 
entering  of  values  for  any  of  these  fields  will  then  set  a  filter  to  select 
only  a  portion  of  the  database  for  the  suranary  report.  For  example,  entering 
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System 


Document  Ref.  Menu 


*  Emissions  Factor  Document  Reference  * 


DataBase  Maintenance 
Summary  Report 
Help  Screen  Review 
Quit  ... 


Air  Quality  Modeling  System 


Emission  Factor  Ref.  DataBase 


Ret/Beg/End/Nxt/Prv/Skp/Mod/Add/Cpy/Del/Lst/Filt/Tally/+/ -/Help/Quit 
Retrieve  a  record  by  its  key  fields 

figure  5.  Emissions  Factor  Document  Reference  Menu  (top)  and 
Database  Maintenance  Option  Menu  (bottom). 
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Figure  6.  Emission  Factor  Reference  Document  Entry  Screen  (top)  and 
Text  Input  Screen  (bottom). 


Filter;  Set/Cancel/Quit?  S 


Figure  7.  Emissions  Data  Summary  Menu  Screen  (top)  and  Filter- 
Setting  Screen  (bottom). 


30 


FI  in  the  AIRCRAFT  IDENTIFICATION  field  and  SOj  in  the  POLLUTANT  field  (bottom 
of  Figure  7)  will  cause  the  filter  to  be  set  to  select  only  SOj  records  for 
aircraft  having  FI  as  the  first  two  characters.  After  the  user  spaces  past 
the  last  open  field  on  this  screen,  the  screen  shown  at  the  top  of  Figure  8  is 
displayed.  If  the  "Go"  option  were  selected,  the  report  shown  at  the  bottom 
of  Figure  8  would  be  produced.  Before  listing  the  report,  the  software  asks 
the  user  if  the  report  is  to  be  sent  to  either  the  screen  or  the  printer. 

Information  from  the  Emission  Factor  Reference  Database  file  can  also 
be  printed  using  procedures  like  those  described  previously  for  the  Emissions 
Database.  However,  for  the  Emission  Factor  Reference  Database  only  a  single 
one-character  field,  the  emission  factor  flag,  can  be  used  for  selecting  or 
filtering  records  to  print.  Therefore,  the  only  possibilities  are  to  print  a 
single  record  from  this  file  or  to  print  all  records.  The  report  selection 
menu  for  this  file  is  accessed  by  selecting  "Emissions  Flag  Ref."  from  the 
DATABASE  MAINTENANCE  menu  (Figure  3)  and  then  "Summary  Report"  from  the 
following  menu  (Figure  4). 

3.  Reindexing  the  Emissions  Database 

Periodically,  after  the  user  has  made  a  number  of  revisions  or 
additions  to  the  database,  it  is  advisable  to  reindex  the  emissions  database. 
Reindexing  updates  internal  database  index  files,  making  database 
manipulations  (such  as  sorting  or  searching)  more  efficient.  Reindexing  is 
performed  by  the  database  management  software  and  takes  very  little  time.  All 
that  is  required  by  the  user  to  perform  this  function  is  to  select  the  "Manual 
Reindex  Function"  item  on  the  DATABASE  MAINTENANCE  menu.  The  software  does 
the  rest  and  very  shortly  returns  control  of  the  system  to  the  user. 

B.  PERFORMING  A  MODEL  RUN 

1.  Selecting  Emissions  Data  for  a  Model  Run 

To  initiate  the  MAILS  program,  load  the  program  and  type  "run". 
Selection  of  the  "Run  MAILS  Model"  from  the  introductory  menu  leads  to  the 
"Main  Menu"  from  which  one  chooses  a  number  of  options  by  using  the 
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1  Air  Quality  Modeling  System 

Emissions  Data  Summary 

AIRCRAFT  IDENTIFICATION 

. . > 

F106 

AIRSPEED  - 

- > 

520  (  MPH  ) 

FUEL  RATE  PER  ENGINE--- 

. --> 

8.64  (IK  Lbs/Hr) 

EMISSION  FACTOR  - 

- > 

1.00  (Lbs/lOOOLbs) 

TYPE  OF  POLLUTANT  - 

. --> 

S02 

NUMBER  OF  ENGINES  - 

--->  :  1 

EMISSION  FLAG  --->  :  0 

EMISSION  RATE  . --- 

- > 

0.00  (Lbs/Hr) 

Set 

-  Filter  :  aircraft='Fr 

and  pollutant 

='S02' 

First  Record  <  AlO  405 

Last  Record  >  T47  100 

Go/</>/Fi 1 t/Beg/End/Nxt/Prv/Skp/Lst/Tal ly/Hel p/Qui t 
Print  or  display  the  report 


MULTIPLE  AIRCRAFT  INSTANTANEOUS  LINE  SOURCE 
EMISSIONS  DATABASE  SUMMARY 


Page  1 


Air 

Craft 

Type 

AVG 

Speed 

(MPH) 

Fuel  Rate 

1  Engine 
(1000  Lb/Hr) 

No. 

of 

Engines 

Pollu 

Type 

Emission 

Factor 

(Lb/ 1000  Lbs) 

Emission 

Rate 

(Lb/Hr) 

F106 

520 

8.64 

1 

S02 

1.00 

8.64 

FlllA 

550 

4.93 

2 

S02 

1.00 

9.86 

FlllD 

550 

5.71 

2 

S02 

1.00 

11.42 

FlllE 

550 

4.93 

2 

S02 

1.00 

9.86 

FlllF 

550 

7.16 

2 

S02 

1.00 

14.32 

F14 

550 

7.40 

2 

S02 

1.00 

14.80 

F15 

550 

5.11 

2 

S02 

1.00 

10.22 

F16 

550 

5.11 

1 

S02 

1.00 

5.11 

F18 

550 

7.50 

2 

S02 

1.00 

15.00 

Press  any 

key  to 

continue. . . 

Figure  8.  Emissions  Data  Summary  Screen  After  Filter  Has  Been  Set 
(top)  and  Resulting  Summary  Report  (bottom). 


32 


up/down  arrow  keys  and  the  "Enter"  key  when  a  choice  is  made.  The  13  Items  on 
the  main  menu  are: 

0  Select  Run  Title 

0  Select  Standard  or  Non-standard  Aircraft  Frequency 
0  Choose  New  Pollutant 
0  Select  Output  File  Name 
0  Select  Nixing  Height 
0  Choose  New  Aircraft 
0  Remove  Aircraft 
0  Edit  Aircraft  Information 
0  Execute  Run 
0  View  Output  File 
0  Retrieve  Parameters 
0  Help  Screen 
0  Quit... 

Each  relevant  menu  Item  Is  selected  by  using  the  arrow  and  "Enter"  keys.  The 
user  enters  all  appropriate  Information  for  a  particular  MAILS  run  and  then 
chooses  the  Item  "Execute  Run"  to  run  the  dispersion  model.  After  the 
completion  of  data  entry  under  any  of  the  menu  Items  from  the  Main  Menu, 
control  Is  returned  to  the  Main  Menu  for  the  next  step.  Current  option 
settings  are  displayed  below  the  Main  Menu  items  which  the  user  should 
monitor.  Corrections  may  be  made  to  any  of  the  entries  controlled  by  the  Main 
Menu  by  simply  choosing  the  appropriate  menu  Item  and  then  editing  the 
relevant  Input  screen.  The  following  comments  relate  to  the  sequential 
selection  of  the  main  menu  Items: 

1.  Run  Title  -  Space  Is  allowed  for  a  maximum  of  70  characters. 
Press  "Enter"  after  selection  to  return  to  main  menu. 

2.  Select  Standard  or  Non-standard  A/C  Frequency  -  The  default 
selection  is  "Standard(s)"  which  Is  displayed  at  the  bottom  of 
the  screen.  To  change,  simply  type  "N";  user  Is  returned  to 
main  menu. 
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3. 


Choose  New  Pollutant  -  Default  selection  is  "CO".  Possible 
other  choices  are  N02,  PART,  S02,  ALL  (plus  HC  and  OTHER  in 
nonstandard  mode).  Model  is  run  automatically  for  all  A  (or  5) 
pollutants  if  "ALL"  is  selected.  If  entry  is  mistyped,  user  is 
prompted  for  correction.  Press  "Enter"  to  return  to  main  menu. 

4.  Select  Output  File  Name  -  Default  selection  is  "P"  (Printer). 
Printer  must  be  ready.  If  "F"  (File)  is  chosen,  output  is 
written  to  user  selected  (name).dat.  (Data  to  be  printed  is 
also  stored  in  the  generic,  temporary  file  "model.dat".  Press 
"Enter"  to  return  to  main  menu. 

5.  Select  Mixing  Height  -  Default  selection  is  "5000"  (ft.).  To 
change  value,  type  in  entire  new  number.  (For  example,  to 
change  from  "5000"  to  "4000",  the  entire  quantity  "4000"  must  be 
typed;  simply  editing  out  the  "5"  by  replacement  with  "4"  will 
result  in  a  final  value  of  4  ft.)  Press  "Enter"  to  return  to 
main  menu. 

6.  Choose  New  Aircraft  -  First  screen  is  for  instructions.  Type 
"Esc"  to  return  to  main  menu  or  type  any  other  key  to  view 
aircraft  list.  Select  one  or  more  aircraft  by  using  arrow  keys 
followed  by  "Ctrl -End".  To  exit  the  aircraft  list,  type  "Esc" 
which  returns  to  an  instruction  screen.  Type  any  key  to  view 
and  edit  aircraft  selection  data.  When  all  Release  Altitudes 
and  Flight  Frequencies  are  entered,  press  "Esc"  to  return  to 
main  menu. 

7.  Execute  Run  -  This  choice  results  in  a  calculation  using  the 
parameters  supplied  by  the  previous  options.  Results  are  output 
either  to  the  printer  or  the  output  file  named  in  4. 

8.  View  Output  File  -  This  choice  provides  for  a  view  of  the 
results  of  an  execution.  Use  either  the  arrow  or  Page  Up,  Page 
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Down  keys  to  peruse  the  file.  Press  "Esc"  to  return  to  the  main 
menu. 

9.  Retrieve  Parameters  -  This  choice  results  in  the  current 

parameters  being  set  to  values  recorded  previously  in  a  file 
(name).txt.  If  used,  the  parameters  in  (naroe.txt)  will 
OVERWRITE  current  settings. 

After  the  completion  of  a  dispersion  model  run,  control  is  again 
returned  to  the  Main  Menu  where  some  or  all  of  the  input  data  may  be  changed 
by  again  choosing  the  relevant  menu  items. 

Dispersion  model  run  time  depends  on  the  type  of  computer  and  on  the 
number  of  aircraft  types  selected,  but  it  is  relatively  short.  For  example, 
on  a  12-MHz,  IBM*-compatible  PC-AT,  a  run  for  three  aircraft  types  would  take 
approximately  a  half  minute  from  the  time  the  user  completes  the  data  entry 
portion  of  the  model  run.  If  the  results  are  to  be  written  to  an  internal 
file,  that  file  is  named  (output  file  name).dat.  If  the  results  are  to  be 
written  to  a  printer,  they  will  be  printed  immediately  if  the  printer  is 
ready.  Whether  or  not  the  printer  is  ready,  the  file  to  be  printed  is  also 
stored  in  model. dit.  If  an  output  file  already  exists,  the  user  will  be  asked 
whether  a  new  run  should  append  the  previous  file  or  overwrite  it.  The 
default  choice  is  to  append,  however  the  user  must  be  cautious  so  that 
previous  results  are  not  mistakenly  overwritten. 

Example  model  outputs  are  shown  in  Figures  12,  14,  and  15. 
Interpretation  of  the  model  output  is  explained  in  the  following  section. 

2.  Applying  Results 

Analysis  of  NAAQS  and/or  PSD  Class  I  increment  consumption  is  the 
intended  application  of  the  MAILS  model.  Because  there  are  three  pollutants 
for  which  PSD  Class  I  increments  have  been  established  (see  Appendix  A),  a 
typical  analysis  for  a  Class  I  area  would  require  three  model  runs  (unless  the 
ALL  option  is  selected):  one  each  for  SO2,  NOg,  and  PART.  If  the  user  has 
selected  the  "standard"  operating  mode  for  the  model  runs,  the  application  of 
the  model  results  is  very  simple.  The  model  computes  and  outputs  the  amount 
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of  each  PSD  Class  I  increment  consumed  bj'  the  route  analyzed  and  also  the 
percentage  of  the  allowable  increment  represented  by  the  impact.  Based  on 
precedent  set  by  EPA  regulations  (40  CFR  51.167),  impacts  less  than  5  percent 
of  an  allowable  increment  were  considered  to  be  insignificant  (Reference  1). 

No  further  analysis  is  required  in  such  cases. 

If  the  predicted  incremental  impact  of  any  route  segment  is  over  5 
percent  of  an  allowable  Class  I  increment,  the  analyst  should  first  verify 
that  the  model  input  values  are  not  overconservative  (flight  frequencies  too 
high  or  altitude  too  low).  If  the  appropriate  model  input  data  have  been 
used,  further  analysis  (a  cumulative  impact  assessment)  may  be  required  in 
order  to  determine  the  total  amount  of  Class  I  increment  consumed  by  the 
combined  impacts  of  the  low-level  route  and  other  types  of  pollutant  sources. 
The  user  should  contact  the  appropriate  local,  state,  or  EPA  air  pollution 
control  agency  for  information  on  existing  PSD  Class  I  increment  consumption 
in  such  cases.  If  other  sources  have  consumed  increment  in  the  PSD  Class  I 
area  of  interest,  in  most  cases  the  appropriate  agency  should  have  some  model 
predictions  that  estimate  the  amounts  of  increment  already  consumed. 

Application  of  the  "nonstandard"  model  run  mode  may  be  necessary  if 
air  quality  standards  change  (for  example,  if  a  24-hour  PSD  Class  I  increment 
were  developed  for  NOj)  or  if  for  any  reason  the  analyst  wants  to  evaluate  the 
impact  of  a  pollutant  for  an  averaging  period  not  represented  in  the 
"standard"  run  mode.  Also,  the  MAILS  model  can  be  used  to  evaluate  impacts  of 
other  pollutants  and  aircraft  types  that  have  not  been  entered  in  the 
permanent  database  for  any  of  the  averaging  periods  (1-hour,  3-hour,  8-hour, 
24-hour,  and  annual)  represented  in  the  "nonstandard"  mode.  This  evaluation 
can  be  accomplished  by  entering  OTHER  for  the  pollutant  type  and  by  entering 
the  appropriate  values  on  the  aircraft  data  entry  screen  (Figure  11). 

C.  EXAMPLE  APPLICATIONS 

The  following  examples  show  the  types  of  data  that  need  to  be  assembled 
before  a  modeling  exercise;  they  illustrate  how  the  model  user  can  focus  the 
air  quality  impact  analysis  on  the  MTR  segments  of  greatest  concern. 


Air  Quality  Modeling  System  MAILS  Model  Main  Menu 


*Air  Quality  Model  Interface:  Main  Menu* 

SELECT  RUN  TITLE 

SELECT  STANDARD/NONSTANDARD  FREQUENCY 

CHOOSE  NEW  POLLUTANT 

SELECT  OUTPUT  FILE  NAME 

SELECT  MIXING  HEIGHT 

CHOOSE  NEW  AIRCRAFT 

REMOVE  AIRCRAFT 

EDIT  AIRCRAFT  INFORMATION 

EXECUTE  RUN 

VIEW  OUTPUT  FILE 

RETRIEVE  PARAMETERS 

HELP  SCREEN 

Quit. . . 

CURRENT  OPTION  SETTING: 

TITLE:  EXAMPLE  1 

FREQUENCY:  S  POLLUTANT:  S02 

OUTPUT  FILE:  EXAMP1.DAT  NO.  AIRCRAFT:  4  MIX  HEIGHT:  5000 


AIRCRAFT 

AlO 

A37 

A4 

A6 

A7 

AV8 

BIB 

B52G 

B52H 

CI30 

C141 

C5 

F106 

FlllA 

FlllD 

FlllE 

FIllF 


BROWSE 


||<C:>||MAILS  IIRec:  1/225 

View  and  Edit  Fields. 


Figure  9.  MAILS  Model  Main  Menu  (top)  and  Aircraft  Selection  Screen  (bottom). 
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Figure  10.  Hypothetical  Nllltary  Training  Route  SegMnts 


1 .  Exampi e  1 


Figure  10  shows  a  hypothetical  set  of  MTR  segments  that  pass  over  a 
PSD  Class  I  area.  Segments  A  and  6  are  portions  of  Individual  routes  that 
join  to  form  a  concurrent  route  segment  (C)  over  the  Class  I  area.  Although 
three  separate  line  sources  exist  over  different  parts  of  the  Class  I  area.  It 
Is  obvious  that  segment  C  would  generate  the  greatest  ground-level  Injects 
because  It  would  have  the  greatest  amount  of  traffic  and  because  the  minimum 
altitudes  of  the  three  route  segments  are  equal.  Therefore,  only  segment  C 
would  be  modeled  to  determine  the  maximum  Impact  with  respect  to  PSD  Class  I 
Increments. 

Table  5  summarizes  the  aircraft  altitude  and  frequency  data  that  would 
be  used  for  estimating  the  air  quality  Impacts  of  segment  C.  The  other 
necessary  model  Input  data  for  these  aircraft  are  selected  from  the  emissions 
database.  This  hypothetical  example  assumes  that  there  are  no  limitations  on 
the  mix  of  the  four  aircraft  types  that  would  be  flown  over  segment  C  In  any 
given  3-hour  or  24-hour  period.  Therefore,  It  Is  important  to  assess  the 
Impacts  of  the  worst-case  aircraft  type  with  regard  to  these  short-term 
periods.  The  B-52H,  the  most  frequently  used  aircraft  on  the  route,  has 
emissions  densities  of  SO2  and  PM-10  (for  which  24-hour  or  3-hour  estimates 
are  necessary)  that  are  roughly  the  same  as,  or  greater  than,  those  of  the 
other  aircraft.  Therefore,  the  maximum  24-hour  and  3-hour  flight  frequencies 
are  based  on  the  assumption  that  all  flights  during  these  periods  could  be 
made  by  6-52H  aircraft.  Because  the  emissions  densities  of  various  aircraft 
can  be  compared  by  browsing  through  the  emissions  database  (see  Figure  9, 
top),  the  3-hour  and  24-hour  frequencies  are  Input  as  zero  for  the  other 
aircraft  types  when  the  model  runs  are  made  for  SOg.  Based  on  scheduling  or 
other  limitations.  It  Is  assumed  that  a  B-52H  could  cross  the  route  segment  a 
maximum  of  6  times  In  3  hours  or  16  times  In  24  hours. 

The  analysis  of  PSD  Class  I  Increment  consumption  Is  accomplished  with 
three  model  runs  (one  if  ALL  Is  selected):  one  each  for  NOj,  SOg,  and  PM-10. 
The  model  Input  and  output  for  the  SO2  run  are  discussed  for  this  example. 

The  MAILS  system  is  started  by  typing  "run"  (displays  title  screen)  or 
"malls"  (bypasses  title  screen)  at  the  DOS  prompt  and  pressing  "Enter."  After 
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TABLE  5.  AIRCRAFT  INPUT  DATA  FOR  EXAMPLE  1 


Segment 

Aircraft 

Minimum 

Altitude  (ft.) 

Maximum  Fliaht  Freouencv 

3-h  24-h  Annual 

A 

R52G 

400 

200 

A 

B52H 

400 

500 

A 

BIB 

400 

100 

B 

FBI  11 

400 

300/450“ 

B 

BIB 

400 

100/150“ 

C 

B52G 

400 

6 

16 

200 

C 

B52H 

400 

500 

C 

BIB 

400 

200 

C 

FBlll 

400 

300 

“The  numerator  denotes  the  total  annual  number  of  flights  (sorties)  on  this 
NTR  segment,  but  because  half  of  the  flights  are  expected  to  make  a  second 
pass  on  segment  B  via  the  "racetrack",  the  total  rt'mber  of  aircraft  passes  for 
segment  B  would  be  as  shown  in  the  denominator. 
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"Run  MAILS  Model"  is  selected  from  the  master  menu,  the  system  displays  the 
model  main  menu  (Figure  9,  top)  to  allow  selection  of  aircraft,  run  mode, 
pollutant  (the  user  selects  SO2),  and  other  data  as  described  in  section 
IV.B.l  above.  When  all  appropriate  input  data  has  been  entered,  the  user 
selects  "Execute  Run"  from  the  main  menu  to  perform  calculations  for  all 
selected  aircraft.  A  confirmation  screen  appears  to  remind  the  user  that 
output  will  be  written  to  either  the  p»*in*er  or  to  the  user-specified  output 
file.  The  user  presses  "enter",  and  model  execution  begins.  The  MAILS  model 
computes  the  maximum  increases  (increments)  in  3-hour,  24-hour  and  annual  SO^ 
concentrations  and  also  the  percentages  of  the  allowable  PSD  Class  I  and  NAAQS 
increments  represented  by  these  impacts.  These  values  are  then  output 
directly  to  the  printer  or  to  a  user-named  print  file,  depending  on  the  user's 
preference  as  selected  during  the  model  input  phase.  Model  output  is  also 
sent  to  a  backup  print  file  (filename  =  MODEL. PRT)  when  the  option  is  chosen 
to  send  output  directly  to  the  printer. 

The  MAILS  results  shown  in  Figure  12  indicate  that  the  maximum  3-hour, 
24-hour  and  annual  SOg  concentrations  are  less  than  5  percent  of  the  allowable 
PSD  Class  I  increments.  Therefore,  the  impacts  of  the  MTR  segment  on  PSD 
Class  I  SO2  increments  are  considered  negligible,  and  no  further  analysis  is 
required  for  SO2.  Additional  model  runs  for  NO2  and  PM-10  would  complete  the 
analysis  of  air  quality  for  the  PSD  Class  I  area. 

2.  Example  2 

Figure  13  shows  a  hypothetical  set  of  MTR  segments,  two  of  which  cross 
two  separate  PSD  Class  I  areas.  Table  6  shows  the  aircraft  altitude  and 
frequency  data  used  for  this  example.  The  other  necessary  model  input  data 
for  these  aircraft  were  those  contained  in  the  emissions  database. 

The  most  efficient  analytical  approach  would  be  to  first  attempt  to 
determine  by  inspection  if  one  of  the  Class  I  area  segments  would  clearly 
cause  a  greater  impact  than  the  other.  The  worst-case  segment  could  then  be 
modeled,  and  if  the  predicted  impacts  were  negligible  (<5  percent  of  any  PSD 
Class  I  increment)  for  all  pollutants,  the  other  segment  would  not  have  to  be 
modeled.  In  this  case,  segment  C  cleariy  has  the  greatest  cumulative 
emissions  density,  while  segment  A  has  a  lower  minimum  altitude.  Therefore, 
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AIRCRAFT 


B526 
B52H 
BIB 
FBI  11 


BROWSE 


SPEED  EMISS_RATE  REL_ALT  FREQJHR  FREQ_24HR  FREQ_ANN  POLLUTANT 


400 

■ 

53.52 

jjj^  1 

6 

16 

200 

S02 

400 

49.92 

Vu  1 

0 

0 

500 

S02 

610 

20.44 

Vn  T 

0 

0 

200 

S02 

550 

■ 

11.42 

JBm 

0 

_ 

0 

300 

S02 

||<C:  >11  MODEL 


IIRec:  1/6 

View  and  Edit  Fields. 


II  II 


Figure  11.  Aircraft  Data  Entry/Edit  Screen. 
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HAILS  -  VERSION  3.0:  MULTIPLE  AIRCRAFT  INSTANTANEOUS  LINE  SOURCE  MODEL 


EXAMPLE  1 


Pollutant  : 
Avg.  Period: 

S02 

3-hour 

No.  of  Aircraft  (Types) 
Mixing  Height 

Aircraft 

Altitude 

(ft) 

Airspeed 

(mph) 

Emiss.  Rate  Flight 
(Ib/hr)  Freq. 

B52G 

400 

400 

53.52  6 

B52H 

400 

400 

49.92  0 

BIB 

400 

610 

20.44  0 

FBlll 

400 

550 

11.42  0 

Total  3-hour  cone. 

****-***** 

4 

5000  ft. 

3 -hour  Cone. 
(micrograffls/m**3) 


.2269 

O.OOE+00 

O.OOE+OO 

O.OOE+00 


.2269 


The  total  3-hour  cone,  is  .9076  %  of  the  PSD 
Class  I  3-hour  increment  for  S02  {  25  micrograms/m**3) 


EXAMPLE  1 

Pollutant  :  S02 

Avg.  Period:  24-hour 


No.  of  Aircraft  (Types) 
Mixing  Height 


Aircraft 

Altitude 

(ft) 

Airspeed 

(mph) 

B52G 

400 

400 

B52H 

400 

400 

BIB 

400 

610 

FBlll 

400 

550 

Emiss.  Rate 
(Ib/hr) 

49.92 

20.44 

11.42 


FI ight 
Freq. 

"'ll' 

0 

0 

0 


ifkititiciticfcic 

4 

5000  ft. 

24-hour  Cone, 
(mi crograms/m**3 ) 

. ’0378 . 

O.OOE+00 

O.OOE+00 

O.OOE-hOO 


Total  24-hour  cone.  = 

The  total  24-hour  cone,  is  .7563  %  of  the  PSD 
Class  I  24 -hour  increment  for  S02  (  5  micrograms/m**3) 


.0378 


EXAMPLE  1 

Pol 1 utant  :  S02 
Avg.  Period:  Annual 


********* 


No.  of  Aircraft  (Types) 
Mixing  Height 


Ai rcraft 

Altitude 

(ft) 

Airspeed 

(mph) 

Emiss.  Rate 
(Ib/hr) 

Flight 

Freq. 

B52G 

400 

400 

53.52 

200 

B52H 

400 

400 

49.92 

500 

BIB 

400 

610 

20.44 

200 

FBlll 

400 

550 

11.42 

300 

4 

5000  ft. 

Annual  Cone, 
(mi crograms/m**3 ) 

.0005 

.0012 

.0001 

.0001 


Total  annual  cone.  = 


.0020 


The  total  annual  cone,  is  .0988  %  of  the  PSD 
Class  I  annual  increment  for  S02  (  2  micrograms/m**3) 


Figure  12.  MAILS  SO2  Results  for  Example  1. 
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TABLE  6.  AIRCRAFT  INPUT  DATA  FOR  EXAMPLE  2 


Hinimuir  Maximum  Flight  Frequency 

Segment  Aircraft  Altitude  (ft.)  3-h  24-h  Annual 


A 

AlO 

200 

400 

A 

F4E 

200 

8 

20 

400 

A 

F16 

200 

400 

B 

F4E 

400 

800 

C 

AlO 

400 

400 

C 

F4E 

400 

12 

30 

1200 

C 

F16 

400 

400 

45 


1t  cannot  readily  be  determined  beforehand  which  segment  would  cause  the 
greatest  ground-level  concentrations  (i.e.,  both  segments  must  be  modeled). 

The  input  data  entry  screens  are  not  shown  for  this  example  since  they 
will  be  basically  the  same  as  for  the  preceding  Example  1  except  for  different 
data  values.  Only  the  output  results  for  S02,  which  are  shown  in  Figures  14 
and  15,  are  discussed  for  this  example. 

The  MAILS  results  indicate  that  the  incremental  3-hour,  24-hour,  and 
annual  SO2  concentrations  will  be  negligible  (<5  percent)  with  respect  to  the 
corresponding  PSD  Class  I  increments  at  both  of  the  hypothetical  PSD  Class  I 
areas  shown  in  Figure  13.  Note  that  the  predicted  concentrations  were 
somewhat  greater  for  segment  A  than  for  segment  C.  Thus,  the  lower  minimum 
altitude  on  segment  A  more  than  outweighed  the  greater  emissions  density  (more 
flights)  on  segment  C.  Additional  model  runs  for  NO2  and  PM-10  for  each  of 
the  Class  1  areas  would  complete  the  air  quality  analysis. 
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MAILS  -  VERSION  3.0:  MULTIPLE  AIRCRAFT  INSTANTANEOUS  LINE  SOURCE  MODEL 


EXAMPLE  2,  SEGMENT  A 


Pollutant 
Avg.  Period: 

Ai rcraf t 

S02 

3-hour 

Altitude 

(ft) 

Ai rspeed 
(mph) 

Mo.  of  Aircraft 
Nixing  Height 

Emiss.  Rate 
(Ib/hr) 

(Types) 

Flight 

Freq. 

AlO 

200 

405 

1.84 

0 

F4E 

200 

550 

14.00 

8 

F16 

200 

550 

5.11 

0 

Total  3-hour 

cone. 

********* 


3 

5000  ft. 

3 -hour  Cone. 
(micrograms/m**3) 

6’6oe+oo 

.1692 

O.OOE+00 


.1692 


The  total  3-hour  cone,  is  .6767  %  of  the  PSD 
Clrss  I  3-hou<  increment  for  S02  (  25  micrograms/ro**3) 


EXAMPLE  2,  SEGMENT  A 


********* 


Pol 1 utant  : 
Avg.  Period: 

Aircraft 


AlO 

F4E 

F16 


S02 

24-hour 

Altitude 

(ft) 

"200 

200 

200 


No.  of  Aircraft  (Types) 
Nixing  Height 


Airspeed 

(mph) 

“405” 

550 

550 


Emiss.  Rate 
(Ib/hr) 

. I’ei" 

14.00 

5.11 


Flight 

Freq. 

..... 

20 

0 


5000  ft. 

24-hour  Cone. 
(micrograms/m**3) 

o!ooe+oo 

.0264 

O.OOE+00 


Total  24-hour  cone.  = 


.0264 


The  total  24-hour  cone,  is  .5286  %  of  the  PSD 
Class  I  24-hour  increment  for  S02  (  5  micrograms/m**3) 


EXAMPLE  2,  SEGMENT  A 


********* 


IWIIWtWIIV  • 

Avg.  Period: 

Aircraft 

Annual 

Altitude 

(ft) 

Airspeed 

(mph) 

Nixing  Height 

Emiss.  Rate 
(Ib/hr) 

-  \  ^  9 - § 

Flight 

Freq. 

:  5000  ft. 

Annual  Cone . 
(micrograms/m**3) 

AlO 

200 

405 

1.84 

400 

.0001 

F4E 

200 

550 

14.00 

400 

.0006 

F16 

200 

550 

5.11 

400 

.0002 

Total  annual  cone.  = 

The  total  annual  cone,  is  .0447  %  of  the  PSD 
Class  I  annual  increment  for  S02  (  2  micrograms/m**3) 


.0009 


Figure  14.  NAILS  SOj  Results  for  Example  2,  Segnent  A. 
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MAILS  - 

VERSION  3.0: 

MULTIPLE  AIRCRAFT  INSTANTANEOUS  LINE 

SOURCE  MODEL 

EXAMPLE  2,  SEGMENT  C 

********* 

Pollutant  : 

S02 

No.  of  Aircraft 

(Types)  : 

3 

Avg.  Period: 

3-hour 

Mixing  Height 

5000  ft. 

Ai rcraf t 

Altitude 

Airspeed 

Effliss.  Rate 

Flight 

3 -hour  Cone. 

(ft) 

(mph) 

(Ib/hr) 

Freq. 

(mi crograms/ra**3 ) 

AlO 

400 

405 

1.84 

0 

O.OOE+00 

F4E 

400 

550 

14.00 

12 

.0863 

F16 

400 

550 

5.11 

0 

O.OOE+00 

Total  3 -hour  cone.  =  .0863 


The  total  3-hour  cone,  is  .3453  %  of  the  PSD 
Class  I  3-hour  Increment  for  S02  {  25  roicrograms/m**3) 


EXAMPLE  2,  SEGMENT  C 


•kitie'kitit’kltic 


Pollutant  : 
Avg.  Period: 


S02 

24-hour 


No.  of  Aircraft  (Types) 
Mixing  Height 


3 

5000  ft. 


Aircraft 

Altitude 

Airspeed 

Eroiss.  Rate 

Flight 

24-hour  Cone. 

(ft) 

(mph) 

(Ib/hr) 

Freq. 

(mi crograms/m**3 ) 

AlO 

400 

405 

1.84 

0 

O.OOE+00 

F4E 

400 

550 

14.00 

30 

.0135 

F16 

400 

550 

5.11 

0 

O.COE+OO 

Total  24-hour  cone.  = 

.0135 

The  total  24 

-hour  cone 

.  is  .2698  %  of  the 

PSD 

Class  I  24-hour 

i ncrement 

for  S02  (  5  micrograms/m**3) 

EXAMPLE  2,  SEGMENT  C 

********* 

Pollutant 

:  S02 

No.  of  Aircraft  (Types) 

:  3 

Avg.  Period 

:  Annual 

Mixing  Height 

:  5000  ft . 

Aircraft 

Altitude 

Airspeed 

Emiss.  Rate 

Flight 

Annual  Cone. 

(ft) 

(mph) 

(Ib/hr) 

Freq. 

(micrograms/m**3) 

AlO 

400 

405 

1.84 

400 

3.52E-05 

F4E 

400 

550 

14.00 

1200 

.0006 

F16 

400 

550 

5.11 

400 

7.19E-05 

Total  annual  cone.  = 


.0007 


The  total  annual  cone,  is 
Class  I  annual  increment  for  S02 


.0349  %  of  the  PSD 
(  2  micrograms/m**3) 


Figure  15.  MAILS  SO2  Results  for  Example  2,  Segment  C. 


48 


SECTION  V 


CONCLUSIONS 


This  user's  guide  provides  a  detailed  description  of  the  interactive  NAILS 
air  quality  model  and  instructions  for  its  use.  MAILS  is  a  siiiq)1e  screening 
model,  providing  estimates  of  worst-case  concentrations  for  NTRs.  The  model 
incorporates  an  aircraft  pollutant  emissions  database  to  allow  convenient  user 
input  for  problem  runs. 

The  intended  application  of  the  MAILS  model  is  the  prediction  of  ground- 
level  pollutant  concentrations  resulting  from  low-flying  (under  3000  feet 
above  ground-level)  military  aircraft  operations  along  a  prescribed  route. 
Previous  studies  have  demonstrated  that  air  quality  impacts  from  low-flying 
military  aircraft  are  (1)  negligible  with  respect  to  NAAQS  and  PSD  Class  II 
air  quality  increments  and  (2)  potentially  significant  with  respect  to  PSD 
Class  I  air  quality  increments,  which  apply  primarily  to  certain  national 
parks  and  wilderness  areas.  Therefore,  the  instructions  and  example 
applications  described  in  this  report  focus  on  the  use  of  the  NAILS  model  for 
analysis  of  air  quality  impacts  on  PSD  Class  I  areas. 
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APPENDIX  A 


NATIONAL  AMBIENT  AIR  QUALITY  STANDARDS  (NAAQS) 

AND  PREVENTION  OF  SIGNIFICANT  DETERIORATION  (PSD)  INCREMENTS* 


Averaging 

PSD  Increments 

Pol 1 utant 

time 

NAAQS 

Class  II 

Class  I 

Nitrogen 

Annual 

100 

25 

2.5 

dioxide 

Sulfur 

3-hour 

1,300^* 

512” 

25” 

dioxide 

24-hour 

365“" 

9l” 

5” 

Annual 

80 

20 

2 

Particulate 

24 -hour 

150** 

30”'” 

8”'” 

matter' 

Annual 

50 

17” 

4” 

Carbon 

1-hour 

40,000” 

monoxide 

8-hour 

10,000” 

Ozone 

1-hour 

235* 

Lead 

Calendar 

quarter 

1.5 

concentrations  are  in  units  of  micrograms/m’. 

‘’Not  to  be  exceeded  more  than  once  per  year. 

'Particulate  matter  under  10  microns  in  diameter  (PM-IO). 
"Listed  PM-10  Class  II  and  I  increments  were  recently  proposed, 
pending  final  rule  promulgation. 

*Not  to  be  exceeded  on  more  than  one  day  per  year. 
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APPENDIX  B 


PREVENTION  OF  SIGNIFICANT  DETERIORATION 
CLASS  I  AIR  QUALITY  AREAS 

Figure  B-1  shows  the  PSD  Class  I  areas  that  were  designated  under  the 
Clean  Air  Act  amendments  of  1977.  A  listing  of  these  areas  by  states  is 
provided  in  Table  B-1.  An  additional  area,  the  Northern  Cheyenne  Indian 
Reservation  in  Montana,  was  redesignated  later  to  PSD  Class  I  status.  More 
detailed  maps  of  particular  PSD  Class  I  areas  can  be  obtained  from  the  U.S. 
Department  of  Interior,  Bureau  of  Land  Management.  The  "state  wilderness 
maps"  available  from  this  agency  should  provide  sufficient  resolution  of  parks 
and  wilderness  areas  for  those  involved  in  the  analysis  of  MTR  air  quality 
impacts  on  PSD  Class  I  areas. 
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TABLE  B-1.  PSO  CLASS  I  AREAS  DESIGNATED  UNDER 
THE  CLEAN  AIR  ACT  ANENDNENTS  OF  1977 


Area  names^ 

State 

Acreage 

Establishing 
Public  Law 

Federal  Land 
Manager 

Sipsey  Wild. 

AL 

12,646 

93-622 

USDA-FS 

Bering  Sea  Wild. 

AK 

41,113 

91-622 

USDI-FWS 

Mount  McKinley  NP 

AK 

1,939,493 

64-353 

USDI-NPS 

Simeonof  Wild. 

AK 

25,141 

94-557 

USDA-FWS 

Tuxedni  Wild. 

AK 

6,402 

91-504 

USDI-FWS 

Chiricahua 

AZ 

9,440 

94-567 

USDI-NPS 

National 

Monument  Wild. 

Chiricahua  Wild 

AZ 

18,000 

88-577 

USDA-FS 

Galiuro  Wild. 

AZ 

52,717 

88-577 

USDA-FS 

Grand  Canyon  NP 

AZ 

1,176,913 

65-277 

USDI-NPS 

Mazatzal  Wild. 

AZ 

205,137 

88-577 

USDA-FS 

Mount  Ba1dy  Wild. 

AZ 

6,975 

91-504 

USDA-FS 

Petrified  Forest 

AZ 

93,493 

85-358 

USDI-NPS 

NP 

Pine  Mtn.  Wild, 

AZ 

20,061 

92-230 

USDA-FS 

Saguaro  Wild. 

AZ 

71,400 

94-567 

USDI-NPS 

Sierra  Ancha  Wild. 

AZ 

20,850 

88-577 

USDA-FS 

Superstition  Wild. 

AZ 

124,117 

88-577 

USDA-FS 

Sycamore  Canyon 

AZ 

47,757 

92-241 

USDA-FS 

Wild. 

Caney  Creek  Wild, 

AR 

14,344 

93-622 

USDA-FS 

Upper  Buffalo 

AR 

9,912 

93-622 

USDA-FS 

Wild. 

Agua  Tibia  Wild. 

CA 

15,934 

93-632 

USDA-FS 

Caribou  Wild. 

CA 

19,080 

88-577 

USDA-FS 

Cucamonga  Wild. 

CA 

9,022 

88-577 

USDA-FS 

Desolation  Wild. 

CA 

63,469 

91-82 

USDA-FS 

Dome  Land  Wild. 

CA 

62,206 

88-577 

USDA-FS 

Emigrant  Wild. 

CA 

104,311 

93-632 

USDA-FS 

Hoover  Wild, 

CA 

47,916 

88-577 

USDA-FS 

John  Muir  Wild. 

CA 

484,673 

88-577 

USDA-FS 

Joshua  Tree  Wild. 

CA 

429,690 

94-567 

USDI-NPS 

Kaiser  Wild. 

CA 

22,500 

94-577 

USDA-FS 

Kings  Canyon  NP 

CA 

459,994 

76-424 

USDI-NPS 

Lassen  Volcanic 

CA 

105,800 

64-184 

USDI-NPS 

NP 

Lava  Beds  Wild. 

CA 

28,640 

92-493 

USDI-NPS 

Marble  Mtn.  Wild. 

CA 

213,743 

88-577 

USDA-FS 
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TABLE  B-1.  PSD  CLASS  I  AREAS  DESIBNATED  UNDER  THE 
CLEAN  AIR  ACT  ANENDNENTS  OF  1977  (CONTINUED) 


Area  names^ 

State 

Acreage 

Establishing 
Public  Law 

Federal  Land 
Manager 

Minarets  Wild. 

CA 

109,484 

88-577 

USDA-FS 

Mokelumne  Wild. 

CA 

50,400 

88-577 

USDA- 

FSPinnacles  Wild. 

CA 

12,952 

94-567 

USDI-NPS 

Point  Reyes  Wild. 

CA 

25,370 

94-544,  94-567USDI-NPS 

Redwood  NP 

CA 

27,792 

90-545 

USDI-NPS 

San  Gabriel  Wild. 

CA 

36,137 

90-318 

USDA-FS 

San  Gorgonio  Wild. 

CA 

34,644 

88-577 

USDA-FS 

San  Jacinto  Wild. 

CA 

20,564 

88-577 

USDA-FS 

San  Rafael  Wild. 

CA 

142,722 

90-271 

USDA-FS 

Sequoia  NP 

CA 

386,642 

26  Stat.  478 
(51st  Cong.) 

USDI-NPS 

South  Warner  Wild. 

CA 

68,507 

88-577 

USDA-FS 

Thousand  Lakes 

CA 

15,695 

88-577 

USDA-FS 

Wild. 

Ventana  Wild. 

CA 

95,152 

91-58 

USDA-FS 

Yolla-Bolly- 

CA 

109,091 

88-577 

USDA-FS 

Middle-Eel  Wild. 

Yosemite  NP 

CA 

759,172 

58-49 

USDI-NPS 

Black  Canyon  of 

CO 

11,180 

94-567 

USDI-NPS 

the  Gunnison 

Wild. 

Eagles  Nest  Wild. 

CO 

133,910 

94-352 

USDA-FS 

Flat  Tops  Wild. 

CO 

235,230 

94-146 

USDA-FS 

Great  Sand  Dunes 

CO 

33,450 

94-567 

USDI-NPS 

Wild. 

La  Garita  Wild. 

CO 

48,486 

88-577 

USDA-FS 

Maroon  Bells- 

CO 

71,060 

88-577 

USDA-FS 

Snowmass  Wild. 

Mesa  Verde  NP 

CO 

51,488 

59-353 

USDI-NPS 

Mt.  Zirkel  Wild. 

CO 

72,472 

88-577 

USDA-FS 

Rawah  Wild. 

CO 

26,674 

88-577 

USDA-FS 

Rocky  Mountain  NP 

CO 

263,138 

63-238 

USDI-NPS 

Weminuche  Wild. 

CO 

400,907 

93-632 

USDA-FS 

West  Elk  Wild. 

CO 

61,412 

88-577 

USDA-FS 

Bradwell  Bay  Wild. 

FL 

23,432 

93-622 

USDA-FS 

Chassahowitzka 

FL 

23,360 

94-557 

USDI-FWS 

Wild. 

Everglades  NP 

FL 

1,397,429 

73-267 

USDI-NPS 

St.  Marks  Wild. 

FL 

17,745 

93-632 

USDI-FWS 

Cohotta  Wild. 

GA 

33,776 

93-622 

USDA-FS 
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TABLE  B>1.  PSD  CLASS  I  AREAS  DES16NATE0  UNDER  THE 
CLEAN  AIR  ACT  ANENDNENTS  OF  1977  (CONTINUED) 


Area  names^ 

State 

Acreage 

Establishing 
Public  Law 

Federal  Land 
Manager 

Okefenokee  Wild. 

GA 

343,850 

93-429 

USDI-FWS 

Wolf  Island  Wild. 

GA 

5,126 

93-632 

USDI-FWS 

Haleakala  NP 

HI 

27,208 

86-744 

USDI-NPS 

Hawaii  Volcanoes 

HI 

217,029 

64-171 

USDI-NPS 

Craters  of  the 

ID 

43,243 

91-504 

USDI-NPS 

Moon  Wild. 

Hells  Canyon 

Wild.^ 

ID 

83,800 

94-199 

USDA-FS 

Sawtooth  Wild. 

ID 

216,383 

92-400 

USDA-FS 

Selway-Jitterroot 

ID 

988,770 

88-577 

USDA-FS 

Wild.^ 

Yellowstone  NP^ 

ID 

31,488 

17  Stat.  32 
(42nd  Cong.) 

USDI-NPS 

Mammoth  Cave  NP 

KY 

51,303 

69-283 

USDI-NPS 

Breton  Wild. 

LA 

5,000-f 

93-632 

USDI-FWS 

Acadia  NP 

ME 

37,503 

65-278 

USDI-NPS 

Moosehorn  Wild. 
(Edmunds  Unit) 
(Baring  Unit) 

ME 

7,501 

(2,782) 

(4,719) 

91-504 

93-632 

USDI-FWS 

Isle  Royal e  NP 

MI 

542,428 

71-835 

USDI-NPS 

Seney  Wild. 

MI 

25,150 

91-504 

USDI-FWS 

Boundary  Waters 

NN 

747,840 

88-577 

USDA-FS 

Canoe  Area  Wild. 

Voyageurs  NP 

NN 

114,964 

99-261 

USDI-NPS 

Hercules-61 ades 

MO 

12,315 

94-557 

USDA-FS 

Wild. 

Mingo  Wild. 

MO 

8,000 

94-557 

USDI-FWS 

Anaconda- Pi ntlar 

MT 

157,803 

88-577 

USDA-FS 

Wild. 

Bob  Marshall  Wild. 

MT 

950,000 

88-577 

USDA-FS 

Cabinet  Mtns. 

MT 

94,272 

88-577 

USDA-FS 

Wild. 

Gates  of  the  Mtn 

MT 

28,562 

88-577 

USDA-FS 

Wild. 

Glacier  NP 

MT 

1,012,599 

61-171 

USDI-NPS 

TABLE  B-1.  PSD  CLASS  1  AREAS  DESIGNATED  UNDER  THE 
CLEAN  AIR  ACT  ANENDNENTS  OF  1977  (CONTINUED) 


Area  names^ 

State 

Acreage 

Establishing 
Public  Law 

Federal  Land 
Manager 

Medicine  Lake 

MT 

11,366 

94-557 

USDI-FWS 

Wild. 

Mission  Mtn.  Wild. 

MT 

73,877 

93-632 

USDA-FS 

Red  Rock  Lakes 

MT 

32,350 

94-557 

USDI-FWS 

Wild. Scapegoat  Wild. 

MT  239, 

295 

92-395 

USDA-FS 

Se 1 way -Bitterroot 

MT 

251,930 

88-577 

USDA-FS 

Wild.^ 

U.L.  Bend  Wild. 

MT 

20,890 

94-557 

USDI-FWS 

Yellowstone  NP^ 

MT 

167,624 

17  Stat.  32 
(42nd  Cong.) 

USDI-NPS 

Jarbridge  Wild. 

NV 

64,667 

88-577 

USDA-FS 

Great  Gulf  Wild. 

NH 

5,552 

88-577 

USDA-FS 

Presidential 

NH 

20,000 

93-622 

USDA-FS 

Range-Dry  River 

Wild. 

Brigantine  Wild. 

NJ 

6,603 

93-632 

USDI-FWS 

Bandelier  Wild. 

NM 

23,267 

94-567 

USDI-NPS 

Bosque  del  Apache 

NM 

30,850 

93-632 

USDI-FWS 

Wild. 

Carlsbad  Caverns 

NM 

46,435 

71-216 

USDI-NPS 

NP 

Gila  Wild. 

NM 

433,690 

88-577 

USDA-FS 

Pecos  Wild. 

NM 

167,416 

88-577 

USDA-FS 

Salt  Creek  Wild. 

NM 

8,500 

91-504 

USDI-FWS 

San  Pedro  Parks 

NM 

41,132 

88-577 

USDA-FS 

Wild. 

Wheeler  Peak  Wild. 

NM 

6,027 

88-577 

USDA-FS 

White  Mtn.  Wild. 

NM 

31,171 

88-577 

USDA-FS 

Great  Smoky  Mtns. 
NP® 

NC 

273,551 

69-268 

USDI-NPS 

Joyce  Kilmer- 

NC 

10,201 

93-622 

USDA-FS 

Slickrock  Wild. 

Linville  Gorge 

NC 

7,575 

88-577 

USDA-FS 

Wild. 

Shining  Rock  Wild. 

NC 

13,350 

88-577 

USDA-FS 

Swanguarter  Wild. 

NC 

9,000 

94-557 

USDI-FWS 

Lostwood  Wild. 

NO 

5,557 

93-632 

USDI-FWS 
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TABLE  B-1.  PSD  CLASS  1  AREAS  DESIGNATED  UNDER  THE 
CLEAN  AIR  ACT  ANENDNENTS  OF  1977  (CONTINUED) 


DvBd  names^ 

State 

Acreage 

Establishing 
Public  Law 

Federal 

Manager 

Theodore  Roosevelt 

NO 

69,675 

80-38 

USDI-NPS 

NMP 

Wichita  Mtns. 

OK 

8,900 

91-504 

USDI-FWS 

Wild. 

Crater  Lake  NP 

OR 

160,290 

57-121 

USDI-NPS 

Diamond  Peak  Wild. 

OR 

36,637 

88-577 

USDA-FS 

Eagle  Cap  Wild. 

OR 

293,476 

88-577 

USDA-FS 

Gearhart  Mtn. 

OR 

18,709 

88-577 

USDA-FS 

Wild. 

Hells  Canyon 

OR 

108,900 

94-199 

USDA-FS 

Wild.*' 

Kalmiopsis  Wild. 

OR 

76,900 

88  577 

USDA-FS 

Mtn.  Lakes  Wild, 

OR 

23,071 

88-577 

USDA-FS 

Mt.  Hood  Wild. 

OR 

14,160 

88-577 

USDA-FS 

Mt.  Jefferson 

OR 

100,208 

90-548 

USDA-rS 

Wild. 

Mt.  Washington 

OR 

46,116 

88-577 

USDA-FS 

Wild. 

Strawberry  Mtn. 

OR 

33,003 

88-577 

USDA-FS 

Wild. 

Three  Sisters 

OR 

199,902 

88-577 

USDA-FS 

Wild. 

Cape  Remain  Wild. 

SC 

28,000 

B-632 

USDI-FWS 

Badlands  Wild. 

SO 

64,250 

94-567 

USDI-NPS 

Wind  Cave  NP 

SO 

28,060 

57-16 

USDI-NPS 

Great  Smoky  Mtns. 
NP^ 

TN 

241,207 

69-268 

USDI-NPS 

Joyce  Kilmer- 

TN 

3,832 

93-622 

USDA-FS 

Slickrock  Wild. 

Big  Bend  NP 

TX 

708,118 

74-157 

USDI-NPS 

Guadalupe  Mtns. 

TX 

76,292 

89-667 

USDI-NPS 

NP 

Arches  NP 

UT 

65,098 

92-155 

USDI-NPS 

Bryce  Canyon  NP 

UT 

35,832 

68-277 

USDI-NPS 

Canyonlands  NP 

UT 

337,570 

88-590 

USDI-NPS 
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TABLE  B-1.  PSD  CLASS  I  AREAS  DESIGNATED  UNDER  THE 
CLEAN  AIR  ACT  ANENDNENTS  OF  1977  (CONTINUED) 


find  names^ 

State 

Acreage 

Establishing 
Public  Law 

Federal 

Manager 

Capitol  Reef  NP 

UT 

221,896 

92-507 

USDI-NPS 

Zion  NP 

UT 

142,462 

68-83 

USOI-NPS 

Lyle  Brook  Wild. 

VT 

12,430 

93-622 

USDAFS 

Virgin  Islands  NP 

VI 

12,295 

84-925 

USDI-NPS 

James  River  Face 

VA 

8,703 

93-622 

USDA-FS 

Wild. 

Shenandoah  NP 

VA 

190,535 

69-268 

USDI-NPS 

Alpine  Lakes  Wild. 

WA 

303,508 

94-357 

USDA-FS 

Glacier  Peak  Wild. 

WA 

464,258 

88-577 

USDA-FS 

Goat  Rocks  Wild. 

WA 

82,680 

88-577 

USDA-FS 

Mount  Adams  Wild. 

WA 

32,356 

88-577 

USDA-FS 

Mount  Rainier  NP 

WA 

235,239 

30  Stat.  993 
{55th  Cong.) 

USDI-NPS 

North  Cascades  NP 

WA 

503,277 

90-554 

USDI-NPS 

Olympic  NP 

WA 

892,578 

75-778 

USDI-NPS 

Pasayten  Wild. 

WA 

505,524 

90-554 

USDA-FS 

Dolly  Sods  Wild. 

WV 

10,215 

93-622 

USDA-FS 

Otter  Creek  Wild. 

WV 

20,000 

93-622 

USDA-FS 

Rainbow  Lake 

WI 

6,388 

93-622 

USDA-FS 

Wild. 

Bridger  Wild. 

WY 

392,160 

88-577 

USDA-FS 

Fitzpatrick  Wild. 

WY 

191,103 

94-567 

USDA-FS 

Grand  Teton  NP 

WY 

305,504 

81-787 

USDI-NPS 

North  Absaroka 

WY 

351,104 

88-577 

USDA-FS 

Wild. 

Teton  Wild. 

WY 

557,311 

88-577 

USDA-FS 

Washakie  Wild. 

WY 

686,584 

92-476 

USDA-FS 

Yellowstone  NP^ 

WY 

2,020,625 

17  Stat.  32 
(42nd  Cong.) 

USDI-NPS 
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TABLE  B-1.  PSD  CLASS  I  AREAS  DESIGNATED  UNDER  THE 
CLEAN  AIR  ACT  ANENDNENTS  OF  1977  (CONTINUED) 


Establishing  Federal 
Brad  names^  State  Acreage  Public  Law  Manager 


Applicable  U.S. 

Area  Name  Province  Acreage  Public  Law 


Roosevelt  Campobello  New  Brunswick  2,721  88-363 

International  Park^  Canada 


^Wilderness  is  abbreviated  as  Wild.,  National  Park  and  NP,  and 
National  Memorial  Park  as  NMP. 

^Hells  Canyon  Wilderness,  193,840  acres  overall,  of  which 
108,900  acres  are  in  Oregon  and  83,800  acres  are  in  Idaho. 

^Selway  Bitterroot  Wilderness,  1,240,618  acres  overall,  of 
which  988,770  acres  are  in  Idaho  and  25,930  acres  are  in  Montana. 

^Yellowstone  National  Park,  2,219,737  acres  overall,  of  which 
2,020,625  acres  are  in  Wyoming,  167,624  acres  are  in  Montana,  and 
31,488  acres  are  in  Idaho. 

^Great  Smoky  Mountains  National  Park,  514,577  acres  overall,  of 
which  273,551  acres  are  in  North  Carolina,  and  241,207  acres  are  in 
Tennessee. 

^Ooyce  Kilmer-Slickrock  Wilderness,  14,033  acres  overall,  of 
which  10,201  acres  are  in  North  Carolina,  and  3,832  acres  are  in 
Tennessee. 

^Section  162(a)  designates  all  international  parks  as  mandatory 
Class  I  areas.  This  designation  indicates  Congressional  intent  to 
prevent  visibility  impairment  from  U.S.  air  pollution  sources. 
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APPENDIX  C 


MAILS  DISPERSION  MODEL  FORTRAN  CODE 


$INCUJDE:  » PORBXBC. INC • 

MULTIPLE  AIRCRAFT  INSTANTANEOUS  LINE  SOURCE  (HAILS)  MODEL 
VERSION:  1.1,  DATED  2/15/90 

AUTHORS:  EDWARD  LIEBSCH,  RESEARCH  ASSOCIATE, 

(615-574-2702) 

R.D.  SHARP,  COMPUTING  SPECIALIST  II 
(615-576-2308) 

OAK  RIDGE  NATIONAL  LABORATORY 

OPERATED  BY  MARTIN  MARIETTA  ENERGY  SYSTEMS,  INC. 
P.O.  BOX  2008 
OAK  RIDGE,  TN.  37831 


CHARACTER* 6  AIRC 

DIMENSION  ISTAB(49) ,WS(49) ,ARR(49) ,SC(6) ,SD(6) 

DIMENSION  AIRC(16) ,AS(16) ,QPRIME(16) ,RHT( 16) , IFR( 5, 16) ,CMAX(16) 

CHARACTER*70  TITLE 

CHARACTER* 4  POLU 

CHARACTER*!  CHAR, STAND, PRDEST 

CHARACTER* 12  FNAME 


DATA  ISTAB  /  7*1,9*2,9*3,14*4,5*5,5*6  / 

DATA  WS  /  0.5, 0.8, 1.0, 1.5, 2. 0,2. 5, 3. 0,0. 5, 0.8, 1.0, 1.5, 2. 0,2. 5, 3.0, 
14.0,5.0,2.0,2.5,3.0,4.0,5.0,7.0,10.0,12.0,15.0,0.5,0.8,1.0,1.5,2.0 
2, 2. 5/3., )./5./7./10., 12 . , IS .,20. ,2. ,2. 5, 3. ,4. ,5. ,2. ,2. 5, 3. ,4. ,5.  / 
DATA  SC,SD  /  24.1667,18.333,12.5,8.333,6.25,4.1667,2.5334,1.8096, 

1  1.0857, .72382, .54287, .36191  / 

OPEN  (UNIT-8, FILE  -  'SAILS. OUT’) 

OPEN (UNIT-1 , FILE- 'MODEL. DAT' ) 

OPEN (UNIT-2 , FI LE- ' TITLE . DAT ' ) 

*************  data  FILE  TITLE.DAT  CONTAINS  A  HEADER  RECORD  (A70)  ** 
READ(2,801)  TITLE 

C  ******************************************************************** 

C 

C  **************  INPUT  DATA  (MODEL.DAT)  ****************************** 
C  AIRC(i)  »  AIRCRAFT  i  IDENTIFIER  (LE  6  CHARACTERS) 

C  AS(1)  -  THE  AIRSPEED  OF  AIRCRAFT  i  (MPH) 

C  QPRIME(i)  -  THE  RELEASE  RATE  OF  POLLUTANT  (LB/HR)  FROM  AIRCRAFT  1 

C  RHT(i)  -  THE  RELEASE  HEIGHT  (FT)  OF  POLLUTANT  FROM  AIRCRAFT  i 

C  AMH  -  THE  SURFACE  MIXING  (INVERSION)  HEIGHT  (FT) 

C  THE  INVERSION  HEIGHT  MUST  BE  GREATER  THAN  THE  RELEASE  HEIGHT 

C  OR  THE  GROUND- LEVEL  CONCENTRATION  IS  ASSUMED  TO  BE  NEGLIGIBLE. 

C  PROGRAM  STOPS  WITH  CODE  1111  IF  A  VALUE  LESS  THAN  OR  EQUAL  TO 

C  THE  RELEASE  HEIGHT  HAS  BEEN  ENTERED. 

C  IFR()c,i)  -  FREQUENCY  OF  FLIGHTS  FOR  AIRCRAFT  i  FOR  TIME  PERIOD  k 
C  k=l(l-HR),  k-2(3-HR),  k-3(8-HR),  k-4 (24-HR) ,  k-5(AN.) 

C  POLU  -  POLLUTANT  NAME:  S02  ,N02  ,PART,  OR  CO  (4  CHARACTERS) 

C  STAND  -  STANDARDS  RUN  FLAG:  S  MEANS  TO  COMPARE  CONCENTRATIONS 

C  AGAINST  STANDARD  PSD  (NAAQS  FOR  1,8-HR)  INCREMENTS 

C  PRDEST  =  FLAG  FOR  PRINT  DESTINATION:  P  MEANS  TO  WRITE  THE  CONC. 

C  TABLE  TO  THE  PRINTER,  AND  TO  FILE  MODEL. PRT 

C  PRDEST  -  F  MEANS  TO  WRITE  THE  CONC. 

C  TABLE  TO  FILE  FNAME  (SEE  NEXT  VARIABLE  DEFINITION) 

C  FNAME  »  FILE  NAME  (LE  12  CHARACTERS)  TO  WHICH  THE  CONC.  TABLE 

C  IS  WRITTEN  IF  PRDEST-F 
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c  ***************************************************************** 

c 

I-O 

1  I-I+l 

READ(1,802,END-15)AIRC(I) , AS (I) , OPRIME{I) ,PMT(I) , AMH, (TFR(K, T) , 
1  K-1 , 5) , POLU, STAND, PRDEST, FNAME 

IF(RHT(I) .GE.AMH)THEN 
WRITE(*,803)RHT(I) ,AMH,I 
STOP  1111 
END  IF 

IF(I.EQ.1)THEN 

AMHFT-AMH 

IF(PRDEST.EQ. •F')THEN 
OPEN (UNIT-9 , FILE- FNAME) 

WRITE(*,'(''  Results  from  MAILS  are  written  to  file:'*, 

1  2X,A12,/) ') FNAME 

ELSE 

OPEN (UNIT-9 , FILE- ' MODEL. PRT • ) 

HRITE(*,'(''  Results  froB  MAILS  are  written  to  file:'', 

1  2X,  "MODEL. PRT"  ,/) ') 

END  IF 
END  IF 
GO  TO  1 
15  NUMAC-I-l 
WRITE (8, 804) 

DO  310  IA-1,NUMAC 
IF(1A.EQ.1)THEN 
WRITE (8, 805)  TITLE 
ELSE 

WRITE(8,806)CHAR(12) , TITLE 
END  IF 

WRITE(8,807)  AIRC(IA) ,RHT(IA) , AMHFT,QPRIME(IA) ,AS(IA) 

WRITE (*, 808) AIRC(IA) 

C  SPECIFY  PORTION  OF  LINE  SOURCE  TO  BE  MODELED  AS  A  DISCRETE  PUFF 
C  IN  UNITS  OF  METERS. 

PL  »  100. 

C  CONVERT  AIRCRAFT  AIRSPEED  FROM  MILES/HR.  TO  METERS/SEC. 

ASM  -  AS (lA) *0.447 

C  CONVERT  RELEASE  HEIGHT  AND  MIXING  HEIGHT  FROM  FEET  TO  METERS 
RHT(IA)  -  RHT(IA) *0.3048 
AMH  -  AMHFT*0.3048 
AMHI  -  1.0/ AMH 

C  CONVERT  EMISSION  RATE  FROM  LB/HR  TO  GRAMS/SEC 
QPRIME(IA)  =  QPRIME(IA) *0.126 
C  COMPUTE  THE  MASS  OF  ONE  PUFF  IN  GRAMS 
Q  -  QPRIME(IA)* PL/ASM 

C  BEGIN  LOOP  OVER  ALL  STABILITY/WIND  SPEED  COMBINATIONS 
DO  200  1-1,49 

C  COMPUTE  THE 'NUMBER  OF  PUFFS  WHICH  PASS  A  CENTERLINE  RECEPTOR 
C  IN  A  ONE-HOUR  PERIOD 

NPUFF  -  WS(I)*3600./PL  +  0.5 
UBARI  =  1.0/WS(I) 

X  -  0.0 
SPSI  -  0.0 
1ST  -  ISTAB(I) 

C  BEGIN  LOOP  OVER  PUFFS 
DO  100  J-1, NPUFF 
PX  =  J  -  0.5 
X  -  PL*PX 
XK  »  .001*X 
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C  COMPUTE  SIGMA-y 

TH  -  .017453293*(SC(IST)-SD(IST)*AL0G(XK)) 

SIGY  »  465.11628*XK*TAN(TH) 

SIGYI  -  l.O/SIGY 

C  CALL  SUBROUTINE  TO  COMPUTE  SZGMA-Z 
CALL  SIGMAZ(XK,SZ,IST) 

SIGZI  >  1.0/SZ 

C  CHECK  IF  SIGMA-Z  IS  LARGE  COMPARED  TO  MIXING  HEIGHT 
C  IF  SO,  ASSUME  UNIFORM  VERTICAL  MIXING 
C 

IF(SZ  .GE.  1.6*AMH)  GOTO  SO 
C  CALCULATE  VERTICAL  TERM  (V) 

V  -  0.0 
A2  «  0.0 

40  VL  -  V 

A2  -  A2  2.0 

HMA2  -  A2*AMH 

A3  -  (HMA2-RHT(IA) )*SIGZI 

A4  -  (HMA2>RHT(IA))*SIGZI 

A3  -  -.5*A3*A3 

A4  -  -.5*A4*A4 

AS  >  0.0 

IF(A3  .GT.  -38.)  A5  -  EXP(Aj) 

A6  •  0.0 

IF(A4  .GT.  -38.)  A6  >  EXP(A4) 

V  -  V  +  AS  +  A6 

IF{ABS(V-VL)  .GT.  I.E-8)  GOTO  40 
A2  -  -.S*RHT(IA)*RHT(IA)*SIGZI*SIGZI 
IF(A2  .GT.  -38.)  V  -  EXP(A2)  +  V 

C  EQUATION  FOR  TOTAL  EXPOSURE  FROM  AN  INSTANTANEOUS  PUFF  RELEASE 
PSI  -  .318309886*Q*SIGyi*SIGZI*UBARI*V 
GOTO  90 

C  EQUATION  FOR  TOTAL  EXPOSURE  FROM  AN  INSTANTANEOUS  PUFF  RELEASE  WITH 
C  UNIFORM  VERTICAL  NIXING 

50  PSI  «  Q*SIGYI*AMHI*UBARI*. 39894228 

SUM  EXPOSURE  FOR  PUFF  AND  GET  ANOTHER  PUFF 
90  SPSI  -  SPSI  +  PSI 
100  CONTINUE 

CONVERT  ONE-HOUR  EXPOSURE  TO  A  ONE-HOUR  CONCENTRATION 
CHI  >  SPSI/3600.0 
ARR(I)  »  CHI 
200  CONTINUE 

WRITE  OUTPUT  TO  FILE 
WRITE (8, 809) 

ISTOLD  ■  1 
DO  300  1-1,49 

IF(ISTAB(I)  '.NE.  ISTOLD)  WRITE(8,810) 

WRITE(8,811)  ISTAB(I) ,WS(I) ,ARR(I) 

CMAX(IA)-AMAX1(CHAX(IA) ,ARR(I)) 

ISTOLD  -  ISTAB(I) 

300  CONTINUE 
310  CONTINUE 

CALL  OUTPUT ( AIRC , AS , QPRIME , RHT , AMHFT , I FR , CMAX , NUMAC , POLU , STAND , 

1  PRDEST, TITLE) 

STOP 
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301  FORHAT(A70) 

802  FORMAT(A6,F3.0,F8.0,F5.0,F5.0,4I3,I4,A4,A1,A1,A12) 

803  FORMATC  The  specified  RELEASE  HEIGHT  ',F6.0,'  is  higher  than  the 

1  MIXING  HEIGHT  ',F7.0,/,'  for  AIRCRAFT  #M3) 

804  FORMAT ( 2 5X, 'MAILS  -  VERSION  1.1  (2/15/90) 25X, 'MULTIPLE  AIRCRAF 
IT  INSTANTANEOUS' ,/, 2 5X, 'LIME  SOURCE  MODEL',///) 

805  FORMAT(32X, '****  TITLE  **** ' ,//,A70,//) 

806  F0RMAT(A1,32X, '****  TITLE  ****',// ,K70 , //) 

807  FORMATC  INPUT  DATA  FOR  AIRCRAFT  '.AS, 

<  /.  - ',14X,' - *,/,*  RELEASE  HEIGHT  -',F10.3 

I, '  FEET',/,'  INVERSION  HEIGHT  -',F10.3,'  FEET',/,'  EMISSION  RATE 

2  -',F10.3,'  LB/HR',/,'  AIRCRAFT  AIRSPEED  -', FIO . 3 , «  MILES/HR'// 
3/) 

808  FORMATC  EXECUTION  CONTINUING .  FOR  AIRCRAFT  ',A6) 

809  FORMAT (IIX, 'STABILITY' ,6X, 'WIND  SPEED' , 5X, 'ONE -HR.  COMC. • ,/, 27X, ' ( 
IM/SEC. ) ’ , 6X, ' (GRAMS/M**3) • ) 

810  FORMATC  ') 

811  F0RMAT(14X,I2,11X,F6.2,8X,E11.4) 

END 

SUBROUTINE  SIGMAZ  CALCULATES  THE  VERTICAL  STANDARD  DEVIATION  OF 
THE  PUFF  CONCENTRATION  DISTRIMITION.  THE  COEFFICIENTS  WERE  TAXEN 
FROM  THE  EPA  ISCST  MODEL  (DATED  88207) . 

SUBROUTINE  SIGMAZ (X, SZ , 1ST) 

DIMENSION  SASIGZ(38) ,SBSIGZ(38) , XI (10, 6) ,INDSGZ(6) 

DATA  SASIGZ  /  122.8, 

1  158.08,170.22,179.52,217.41,258.89,346.75,2*453.85, 

2  90.673,98.483,109.3,61.141,34.459,32.093,32.093,33.504,36.65, 

3  44.053,24.26, 

4  23.331,21.628,21.628,22.534,24.703,26.97,35.42,47.61, 

5  15.209,14.457,13.953,13.953,14.823,16.187,17.836,22.651,27.074, 

6  34.219  / 

DATA  SBSIGZ  /.9447, 

X  1.0542,1.0932,1.1262,1.2644,1.4094,1.7283,2*2.1166, 

1  .93198, .98332,1.0971, .91465, .86974, . 81066, . 64403 ,. 60486, . 56589 , 

X  .51179, .8366, 

2  .81956, .75660, .63077, .57154, .50527, .46713, .37615, .29592, 

3  .81558, .78407, .68465, .63227, .54503, . 46490, . 41507 ,. 32681 ,. 27436 , 

4  .21716  / 

DATA  INDSGZ  /O, 9 , 12 , 13 , 19, 28/ 

DATA  XI  /.1,.15,.2,.25,.3,.4,.5,3.11,1.E20,0.,  . 2 , . 4 , 1 . E20 , 7*0 . , 

II. E20,9*0.,  .3,1.,3.,10.,30.,1.E20,4*0.,  . 1 , . 3 , 1 . , 2 . , 4 . , 10 . , 

2  20, ,40. ,1.E20,0. ,  . 2 , .7 , 1. , 2 . , 3 . , 7 . , 15. , 30 . , 60 , , 1 . E20/ 

1  =  1 

IF(IST  .EQ.  3)  GOTO  20 
10  IF(X-X1(1,IST)  ,LE.  0.0)  GOTO  20 
1  =  1  +  1 
GOTO  10 

20  INDXl  =  INDSdz(IST)  +  I 

SZ  =  SASIGZ (INDXl) *X**SBSIGZ (INDXl) 

SZ  =  AMIN1(SZ,5000.) 

RETURN 

END 

SUBROUTINE  OUTPUT ( AIRC , AS , QPRIME , RHT , AMH , I FR , CMAX , NUMAC , POUJ , 

1  STAND, PRDEST, TITLE) 

INTEGER*2  SYSTEM 
CHARACTER*70  TITLE 
CHARACTER  ANS 
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CHARACTER*6  AIRC 
CHARACTER*?  AVT(5) 

CHARACTER  STAND, CHAR, FREEST 
CHARACTER* 4  FOLD 

DIMENSION  AIRC(16) ,AS(16) ,QFRIME(16) ,RHT(16) , IFR(5, 16) ,CMAX(16) 
DIMENSION  FSO(3,4) ,KNTA(4) ,ITM(3,4) ,DIV(5) ,SFACT(5) ,NTF(5) 

DATA  AVT/ ' 1-hour  ' , •3-hour  * , 'S-hour  ' , '24-hour' , 'Annual  '/ 
FSD(IT,IF)  -  FSD  CLASS  1  INCREMENTS  (BlcrograBs/B**3) 

IF>1  FOR  NO,IF-2  FOR  S02,  IF-3  FOR  FART,  IF-4  FOR  CO 

IT-1  FOR  1-HR,  IT-2  FOR  3-HR,  IT-3  FOR  8-HR,  IT-4  FOR  24-HR 
AND  IT-5  FOR  Annual 

NOTE:  FOR  CO,  FSD (1,4)-  NAAQS  1-HR  STANDARD, 

FSD (2, 4)-  NAAQS  8-HR  STANDARD 
-999.  MEANS  NO  VALUE  ASSIGNED  (NA) 

DATA  FSD/2.5,-999. ,-999. ,  25. ,5. ,2.,  8. ,4. ,-999., 

1  40000. , 10000. ,-999./ 

DATA  KNTA/1,3,2,2/ 

DATA  ITM/5,0,0,  2,4,5,  4,5,0,  1,3,0/ 

DATA  NTF/5*0/ 

DATA  SFACT/1.0,0.5,0.33,0.25,0.1/ 

DATA  DIV/1 . 0, 3 .0,8.0, 24.0,8760.0/ 

HRITE(9,800) 

IF(FOLU.EQ. 'N02  ' ) IF-1 
IF ( FOLD. EQ. '802  ■ ) IF-2 
IF ( FOUI . EQ . • FART • ) IF-3 
IF(FOU;.EQ. 'CO  ')IF-4 
IAMH-IFIX(AMH-t-.5) 

IF(STAND.EQ. 'S')THEN 
NUMT-KNTA(IF) 

NLTAB-15 
ELSE 
NUMT-5 
NLTAB-12 
END  IF 

DO  5  NTAB-1,NUMT 

IF ( NTAB . EQ . 1 ) CALL  TABLEC (NUMT , NLTAB , NUMAC , NTF) 

IF(NTF(NTAB) . EQ. 1 .AND. NTAB. NE. 1)THEN 
WRITE(9,8001)CHAR(12) , TITLE 
ELSE 

WRITE (9, 801) TITLE 
END  IF 

WRITE (9,802) FOLU , NUMAC 
IF(STAND.EQ. 'S')THEN 
K-ITM(NTAB,IF) 

ELSE 
K-NTAB 
END  IF 

WRITE ( 9 , 803 ) AVT ( K) , lAMH 
WRITE ( 9 , 804 ) AVT (K) 

HRITE(9,805) 

WRITE (9, 806) 

CT-0.0 

DO  10  IA-1, NUMAC 
C  TO  FRINT  FARAMETERS  AS  INTEGERS 

IRHT-IFIX (RHT ( lA) / . 3048+ . 5) 
lAS-IFIX ( AS ( lA) + . 5 ) 

QF-QFRIME ( lA) / . 126 

C  CONCENTRATION  UNIT  OF  C  IS  MICROGRAMS 

C=( (CMAX(IA) *FLOAT(IFR(K,IA) ) )/DIV(K) ) *SFACT(K) *1.0E+6 
CT-CT+C 
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IF(C.GE.1.0E-4)THEN 

WRITE(9,807)AIRC(IA) ,IRHT,IAS,QP,IFR(K,IA) ,C 
ELSE 

WRITE(9,e08)AIRC(IA) ,IRHT,IAS,QP,IFR(K,IA)  ,C 
END  IF 

10  CONTINUE 

WRITE(9,809) 

IP(CT.GE.1.0E-4)THEN 

IF(AVT(K) .NE. • Annual') WRITE(9, 810) AVT(K) ,CT 
IF(AVT(K) .EQ. 'Annual •)WHITE(9, 811) CT 
ELSE 

IF(AVT(K) .NE. 'Annual')WRITE(9,812)AVT(K) ,CT 
1F(AVT(K) .EQ. 'Annual')WRITE(9,813)CT 
END  IF 

IF(STAND.EQ. 'S')THEN 
PERC-CT/ PS D ( NTAB , I P ) * 10 0 . 

1PERC-IFIX(PSD(NTAB,IP) ) 

1F(PERC.GE. 0.01) THEN 

IF{AVT(K) .NE. 'Annual')WRITE(9,814)AVT(K) ,PERC,AVT(K) .POLD,IPERC 
IF(AVT(K) .EQ. 'Annual')WRITE(9,816)PERC.POLU,IPERC 
ELSE 

IF(AVT(K) .NE. ’Annual')WRITE(9,815)AVT(K) ,PERC,AVT(K) ,POLU,IPERC 
IF(AVT(K) .EQ. 'Annual')WRITE(9,817)PERC,POLU,IPERC 
END  IF 
END  IF 
5  CONTINUE 

800  FORMAT( 'MAILS  -  VERSION  1.1  (2/15/90);  MULTIPLE  AIRCRAFT  INSTANTAN 
lEOUS  LINE  SOURCE  MODEL') 

801  FORMAT (//, A7 0, '  *•*******',/) 

802  FORMAT( IX, 'Pollutant  :  ' ,A4,T30, 'No.  of  Aircraft  (Typaa)  :  ', 

1  15) 

803  FORMAT ( IX, 'Avg.  Period;  ' ,A7,T30, 'Mixing  Height  ;  ', 

1  15 , '  ft . ' , / ) 

804  FORMATCAlrcraft' ,3X, 'Altitude', 3X, 'Airspeed' ,2X, 'Eaiss.  Rate',2X, 

1  'Flight' ,4X,A7, '  Cone.') 

805  FORMAT(T13, ' (ft) ' ,T24, ' (mph) ' ,T35, ' (Ib/hr) ' ,T46, 'Freq. ' ,T59, 

1  'Cone. '  ,T54, '  (]BlerograBS/n**3) ') 

806  FORMATC - ',T12,  ' - ’,T23,  ' - ',T33,  ' - 

807  ’f0RMAT(2X,A6!i8,3X,I8,6X,0PF8.2,5X,I4,7X,0PF8.4) 

808  FORMAT ( 2X , A6 , 18 , 3X, 18 , 6X, 0PF8 . 2 , 5X, 14 , 7X , 1PE8 . 2 ) 

809  FORMAT  (T54,' - ') 

810  FORMAT (T3 3, 'Total  ',A7,'  cone.  =  ',0PF8.4) 

811  FORMAT(T33, 'Total  annual  cone.  «  ',0PF8.4) 

812  FORMAT (T3 3, 'Total  ',A7,'  cone.  =  ',1PE8.2) 

813  FORMAT (T33 , 'Total  annual  cone.  »  ',1PE8.2) 

814  FORMAT{/,9X, 'The  total  ',A7,'  cone,  is  ',F8.4,'  %','  of  the  PSD', 

<  /f 

1  4X, 'Clase  I  ',A7,'  increment  for  ' ,A4, ' ( ' ,15, '  micrograms/B**!) ' ) 

815  FORMAT (/,4x, 'The  total  ',A7,'  cone,  is  ',1PE8.2,'  %','  of  the  PSD 

1  4X, 'Class  I  ',A7,'  increment  for  ' ,A4, ' ( ' ,15, '  micrograms/B**3) ' ) 

816  FORMAT (/,9X, 'The  total  annual  cone,  is  ',F8.4,'  %','  of  the  PSD', 

<  / , 

1  4X, 'Class  I  annual  increment  for  ' ,A4, ' ( ' , 15, '  micrograBS/B**3) ' ) 

817  FORMAT(/,9X, 'The  total  annual  cone,  is  ',1PE8.2,'  %','  of  the  PSD 

1  4X, 'Class  I  annual  increment  for  '  ,A4, '  ( ' ,  15,  '  Bicrograiiis/m**3)  ' ) 
8001  FORMAT ( Al, A7 0, '  ***•*****•,/) 

CLOSE(UNIT=9) 
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c 


c 


*****************  PRIMT  TH£  RESULT  FILE  ?  * 

IF(PRDEST.EQ. 'PMTHEM  _ 

WRITE '  la  the  PRINTER  (LPTl)  ready; 
READ(*, • (Al) •)ANS 
IF(AMS.EQ. *T* .OR.AHS.EQ. *y‘)THEN 
1-SVSTEM( ‘PRINT  MODEL. PRT'C) 


Y  or  N?  "  ) ' ) 


ELSE 

T*SYS*rEM(  *  CIjS  *  C) 

WRITE(*, •(//,'•  Result  file  MODEL. PRT  cannot  bo  printed  by  this 
lprogra»“,/.14X,  "until  LPTl  is  ready!  ")  ') 

WRITE(*,'('*  Hit  ENTER  to  continue.  •*)') 

READ(*, • (Al) ')ANS 
END  IF 


RETURN 

END 

SUBROUTINE  TABLEC(NUMT,NLTAB,NUMAC,NTP) 
DIMENSION  NTP{5) 

NTP(l)  IS  INDICATOR  THAT  TABLE  i  BEGINS 
NLPAGE>60 
IT-0 
NLPB-0 
1  IT-IT+l 

IF(1T.GT.NUMT)G0  TO  10 
NLP-NLTAB+NUMAC+ ( 6-IT) /S+NLPB 
IF (NLP. LE . NLPAGE) THEN 
NLPB-NLP 
GO  TO  I 
ELSE 

NTP(IT)-1 
NLPB-NLTAB+HUMAC 
GO  TO  I 
END  IF 
10  CONTINUE 
RETURN 
END 
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APPENDIX  D 


THE  MULTIPLE  AIRCRAFT  INSTANTANEOUS 
LINE  SOURCE  MODEL  (MAILS) 


BACKGROUND 

Low-altitude  flight  operations  are  a  critical  ingredient  in  maintaining  a 
well -trained,  combat  ready  Air  Force.  However,  operations  along  military 
training  routes  (MTRs),  due  to  their  proximity  to  the  surface  and  repetitive 
flights  within  a  relatively  narrow  corridor,  often  raise  community  concerns 
regarding  adverse  effects  on  visibility,  noise  levels,  and  air  pollutant 
concentrations.  Such  concerns  may  result  in  restrictions  on  acquisition  and 
use  of  airspace,  with  a  negative  impact  on  effective  training  operations. 
Thorough,  defensible  environmental  impact  documentation  addresses  such 
concerns,  and  provides  for  compliance  with  legislation  such  as  the  National 
Environmental  Policy  Act  (NEPA).  Guidance  and  direction  for  preparing 
environmental  documentation  regarding  Air  Force  operations  is  provide  by  AFR 
19-2,  USAF  Environmental  Impact  Analysis  Process  (EIAP).  The  purpose  of  the 
MAILS  model  is  to  support  the  EIAP  regarding  air  quality  impact  of  operations 
along  MTRs.  It  was  developed  to  fill  the  gap  between  simple  lookup  tables  from 
which  total  emissions  along  an  MTR  could  be  estimated,  and  complex  dispersion 
models  requiring  considerable  time  and  effort  in  terms  of  data  gathering, 
input,  and  post-processing  for  applicability  to  MTR  scenarios. 

MAILS 

The  MAILS  model  itself  is  an  interactive  air  quality  screening  tool  that 
combines  a  military  aircraft  engine  emissions  database  with  a  gaussian  puff 
dispersion  model.  MAILS  operates  on  IBM-compatible  PC-AT,  386,  or  486 
microcomputers  with  DOS  version  3.0  or  later.  The  model  requires  360  Kb 
available  RAM,  can  take  advantage  of  a  math  coprocessor  if  one  is  available, 
and  should  normally  be  installed  on  and  run  from  a  hard  disk. 

AIRCRAFT  EMISSIONS  DATABASE 

The  MAILS  emissions  database  contains  information  for  roughly  50  military 
aircraft  types.  For  each  aircraft,  the  database  includes  the  following 
information:  airspeed,  fuel  rate  (per  engine),  emission  factor  (pounds 
pollutant  per  pound  fuel  burned),  number  of  engines,  emissions  flag  (a 
reference  describing  where  emissions  data  for  the  aircraft  was  obtained),  and 
total  emission  rate  (fuel  rate  x  emissions  factor  x  number  of  engines). 
Emissions  data  is  provided  for  five  pollutants  (sulfur  dioxide,  nitrogen 
dioxide,  particulate  matter,  carbon  monoxide,  and  hydrocarbons) ,  and  is  easily 
updated  by  the  user  should  new  information  become  available  or  for  application 
to  special  cases. 
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DISPERSION  MODEL 


The  dispersion  modeling  portion  of  MAILS  utilizes  a  gaussian  puff  dispersion 
model,  approximating  a  line  source  by  dividing  the  emitted  plume  into  puffs, 
each  with  pollutant  mass  equal  to  that  emitted  during  a  100  meter  segment  of 
the  flight  path.  Ground  level  concentrations  at  plume  centerline  are  based  on 
summing  the  contributions  of  all  puffs  that  pass  over  a  fixed  receptor  during 
a  one-hour  time  period.  Consistent  with  the  screening  approach  taken  by  MAILS, 
the  ambient  wind  is  assumed  parallel  to  the  flight  track  and  homogeneous  in 
space  and  time.  In  addition,  a  matrix  of  49  wind  speed/atmospheric  stability 
conditions  are  used  to  calculate  possible  concentrations  at  the  receptor,  from 
which  the  highest  is  selected  for  reporting.  Note  that  the  same  49  atmospheric 
conditions  are  those  used  in  the  EPA-approved  screening  model  PTPLU,  and  that 
the  dispersion  coefficients  are  equal  to  those  used  in  the  EPA  ISCST  model 
("rural"  mode  coefficients  are  used,  since  MTRs  rarely  pass  over  urban  areas). 
Multiple  aircraft  passes  within  one  hour  are  simply  summed  together.  Multiple 
passes  over  longer  time  periods  are  summed  and  then  adjusted  according  to  the 
requested  averaging  time.  MAILS  may  be  used  to  estimate  1-hour,  3-hour,  8- 
hour,  and  24-hour,  and  annual  concentration  averages.  Adjustment  factors  for 
periods  greater  than  one  hour  account  for  variations  in  meteorological 
conditions  and  in  flight  track  positions  relative  to  the  stated  flight 
centerline  location,  and  for  the  intermittent  nature  of  operations  along  MTRs. 

INPUTS/OUTPUTS 

The  user  is  required  to  input  the  aircraft  type(s),  pollutant(s)  of  interest, 
flight  altitude(s),  frequency  of  overflights  of  a  given  ground  location  during 
the  period  of  interest,  and  the  height  of  the  surface-based  temperature 
inversion  (if  applicable).  The  user  is  also  given  default  values  of  airspeed, 
and  emission  rate,  which  he  may  be  adjusted  to  fit  any  special  circumstances. 
The  model  provides  output  (to  printer  or  file)  in  the  form  of  a  table 
detailing  the  component  of  total  concentration  due  to  each  aircraft  type  on  a 
given  route  and  the  total  concentration  of  each  pollutant  type,  both  as  an 
absolute  value  (micrograms/cubic  meter)  and  as  a  percentage  of  NAAQS  and  PSD 
Class  I  standards.  Generally,  if  MTR  emissions  result  in  ground  level 
concentrations  less  than  5  percent  of  the  PSD  Class  I  limits,  the  screening 
analysis  may  be  considered  to  demonstrate  negligible  impact  and  no  further 
analysis  is  necessary.  This  is  usually  the  case  for  MTR  operations. 

VERIFICATION 

Because  data  sets  addressing  the  unique  characteristics  of  MTR  ground 
concentrations  as  modeled  in  MAILS  do  not  exist,  MAILS'  performance  has  been 
tested  by  comparison  with  an  EPA  guideline  model.  The  Industrial  Source 

Complex  -  Short  Term  (ISCST)  model  and  the  MAILS  dispersion  code  were  run  for 

emissions  at  the  same  height  and  rate,  and  under  a  varied  set  of 

meteorological  conditions.  The  results  differed  by  less  than  5  percent  in  all 
cases.  Note  that  configuring  ISCST  for  application  to  an  elevated,  moving  line 
source,  and  inputting  required  aircraft  engine  emissions  data  is  a  laborious 


and  time  consuming  process,  one  of  the  primary  reasons  for  MAILS' 
development.  Details  of  the  model  comparison  are  available  in  the  MAILS  user's 
guide  (ESL-TR-89-59). 

ADVANTAGES  OF  MAILS 

MAILS  offers  several  advantages  over  previous  methods  of  screening  for  air 
quality  impact  along  training  routes.  It  combines  emissions  data,  aircraft 
configuration  data,  and  dispersion  code  into  one  user-friendly  software 
package.  It  is  optimized  to  provide  ground  track  concentration  estimates  below 
training  routes,  where  previous  models  required  laborious  adaptation  to  this 
purpose.  The  MAILS  database  is  readily  updated  as  new  aircraft  types  and 
engine  models  enter  the  operational  inventory.  The  output  from  MAILS  is  given 
in  terms  which  allow  quick  comparison  to  regulatory  standards.  Where 
assumptions  have  been  made  in  the  dispersion  model,  these  have  generally  been 
of  a  conservative  nature,  such  that  the  model  is  unlikely  to  underpredict 
pollutant  concentrations.  An  indication  of  negligible  impact  by  the  model  thus 
has  a  high  degree  of  confidence. 
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